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CHAPTER  I 
INTRODUCTION 

The  disposal  of  dredged  material  has  recently  received  much  atten¬ 
tion.  Section  404(b)  of  the  Federal  Water  Pollution  Control  Act 
Ammendments  of  1972,  P.L.  92-500,  prohibits  discharges  of  dredged 
material  to  navigable  waters  of  the  United  States  unless  permits  are 
issued  through  the  U.S.  Army  Corps  of  Engineers.  In  1975,  guidelines 
on  the  issuance  of  permits  were  published  in  the  Federal  Register. 

Among  the  ecological  impacts  from  dredged  material  disposal  to  be 
examined  were  impairment  of  the  water  column  and  the  covering  of 
benthic  communities.  The  need  for  mathematical  models  to  predict  the 
disposition  of  suspended  solids  resulting  from  disposal  of  dredged 
material  therefore  becomes  apparent. 

The  Corps  of  Engineers  currently  dredges  portions  of  the  upper 
Mississippi  River  to  maintain  a  nine  foot  deep  channel  for  barge  traf¬ 
fic.  The  hydraul ical ly  dredged  material  is  discharged  onto  a  nearby 
island  or  bank  and  the  excess  water  flows  back  into  the  river.  This 
water  contains  suspended  solids,  either  from  the  dredged  sediment  or 
from  the  disposal  site,  and  forms  a  suspended  solids  plume  where  it 
enters  and  rejoins  the  river. 

Much  of  the  modeling  on  suspended  solids  plumes  resulting  from 
dredge  disposal  has  been  for  open  sea  or  estuarine  operations.  Little 
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work  has  been  done  on  dredge  disposal  in  the  river  environment.  The 
"objectives  of  this  study  were  to: 

'  1)  Collect  field  data  on  suspended  solids  and  turbidity  during 
two  dredge  disposal  operations  on  the  Mississippi  River, 

2)  Check  the  utility  of  the  Schubel  and  Carter  (1978)  model  for 
adequately  describing  the  observed  field  data  and  modify,  if 
possible,  to  reflect  river  conditions, 

(.  3)  Examine  other  models  available  to  describe  the  observed  field 
data,  including  the  numerical,  computer  solution  of  Weschler 
and  Cogley  (1977)  (such  models  can  be  used  to  rapidly  generate 
a  number  of  simulations  covering  a  spectrum  of  conditions  ex¬ 
pected  in  the  Mississippi  River),  and 
(  4)  Develop  a  convenient,  analytical  solution  for  the  prediction 
of  suspended  solids  concentrations  caused  by  hydraulically 
dredged  sediment  and  compare  the  model  results  to  field 
measurements. 

The  scope  of  this  modeling  effort  includes  the  utilization  of  ex¬ 
isting  dredge  disposal  mathematical  models,  both  analytical  and 
numerical,  as  well  as  the  development  of  a  new  model.  The  new  model 
is  specifically  derived  for  continuous  nonpoint  source,  sidebank  dis¬ 
posal  type  of  operations  such  as  commonly  practiced  in  the  upper 
Mississippi  River.  Suspended  solids  concentrations  are  predicted. 

This  research  grew  out  of  a  larger  dredging  study  by  a  multi- 
departmental,  multi-disciplinary  consortium  called  the  Great  River 
Environmental  Action  Team,  GREAT  II.  The  GREAT  II  study  reach  of  the 
Mississippi  River  stretches  from  Guttenberg,  Iowa  to  Saverton,  Missouri. 
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CHAPTER  II 
LITERATURE  REVIEW 

Review  of  Models 

Models  for  predicting  the  distribution  of  suspended  solids  result¬ 
ing  from  disposal  of  dredged  material  have  been  proposed  by  Schubel , 


Carter  et  al.,  (1978)  and  Wechsler  and  Cogley  (1977). 

Both  models  be- 

gin  with  the  Fickian  diffusion  equation: 

9C  +  9C_  _  J_  L 

3t  ui  3xi  9x.  ( 

3C 
i  3x.j 

) 

(2.1) 

Rate  of  change 

Rate  of  change  of 

Rate  of  change 

of  suspended  solids 

+ 

suspended  solids  = 

of  suspended 

concentration 

concentration  due 

solids  concen¬ 

to  convection 

tration  due  to 

diffusion 

where  C  refers  to  the  concentration  (mass  per  unit  volume)  of  suspended 
sediment;  refers  to  the  fluid  velocity  in  a  rectangular  coordinate 
system,  x^;  and  K.  refers  to  the  eddy  diffusion  coefficient  in  the  i'th 
direction.  The  models  begin  to  differ  at  this  point  in  the  assumptions 
that  are  made. 


Schubel  and  Carter  Model 

The  model  developed  by  Schubel,  Carter  et  al.,  (1978)  is  for 
estuarine  or  shallow  coastal  dredge  disposal  operations.  The  initial 
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assumptions  are:  1)  the  individual  concentrations  of  the  various  size 

fractions  of  suspended  sediment,  c. ,  can  be  described  by  a  vertically 

D 

averaged  suspended  solids  concentration,  C  =  ^  /  Sc.  dz,  where  D 

o 

is  the  depth  of  the  water  column  and  z  is  vertical  distance  in  the 

Cartesian  coordinate  system,  x,  y,  z;  2)  the  eddy  diffusivities  in  the 

x  and  y  directions,  K  and  K  ,  are  equal  and  independent  of  depth; 

x  y 

3)  the  fluid  velocity  in  the  x  and  y  directions,  u  and  u  ,  are  depth 

x  y 

independent;  and  4)  the  terms  for  vertical  diffusion  and  convection  can 

WC 

be  combined  into  one  term,  -  p-,  where  W  is  the  mean  settling  velocity 


of  the  particles. 


J  £ w.c.dz 
o  11 
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and  wi  is  the  settling  velocity  of  the  individual  particle,  c^.  This 
fourth  assumption  is  based  on  the  assumption  that  the  suspended  solids 
transport  due  to  vertical  diffusion  and  vertical  fluid  velocity  cur¬ 
rents  is  much  smaller  than  the  transport  due  to  the  settling  velocity 
of  the  suspended  solids.  The  resulting  equation  is: 


3C  _  9 _  p  L  „  r  +  L  k  +  L  if  wc 

at  ■  "  9x  xT  ■  9y  v  Sx  Nx  9x  9y  x  9x  '  D 


(2.2) 


Okubo  and  Pritchard  (Okubo,  1962)  proposed  the  solution  assuming 
an  instantaneous  vertical  line  source.  This  solution  is  then  inte¬ 
grated  over  time  to  describe  a  continuous  vertical  line  source.  The 
resulting  equation  is: 


5 


C(x,y,t)  = 


• 

_  — 

2 

^  - 

q 

1 

X  -  uxf 

exp 

.  _X_ 

<i)t ' 

TT  0)2D 

£12  eXP 

cot' 

0 

-  J 

exp 


N 


dt* 


(2.3) 


where  q  is  the  rate  of  suspended  material  added  to  the  plume  (mass  per 

time)  and  co  is  the  diffusion  velocity  (cm/sec).  The  diffusion  velocity, 

a),  is  related  to  the  horizontal  eddy  diffusion  coefficient  by 

Kx  =  co2t.  The  first  and  second  exponential  terms  in  the  integral  refer 

to  diffusion  of  suspended  solids  in  the  x  and  y  direction,  while  the 

third  exponential  term  in  the  integral  represents  particle  settling. 

The  model  is  not  used  in  this  form,  however.  First,  x,  y,  and  t‘ 

are  nondimensionalized  to  x  ,  y  and  t  ,  where 
★ 

X  =  X-  Uxt 

y  -  y*  uxt 

★ 

t1  =  t  t 


The  resulting  equation  is: 

1 


C(x,y,t)  =  — | 


TT«2Dt 


* 


exp 


•  M 


dt 


Tt*P 


exp  - 


★  * 
x  -t 


exp  - 


■  * 
t* 

(2.4) 


The  integral  term  is  defined  as  a  function,  G,  of  x  ,  y*,  w/u  and  y, 
where  y  =  Wt/D  and  relates  the  plume  age,  t,  to  the  settling  time, 

W/D.  Normalizing  Equation  2-4  by  the  concentration  at  the  plume  front 


(at  distance  uxt),  the  final  form  of  the  model  is  obtained: 


★ 

For  the  centerline,  y  =  0,  Equation  2.5  reduces  to: 

The  solution  to  the  model  for  the  plume  centerline  is  contained 

in  a  series  of  graphs  of  V  vs-  x*  with  w/uv  and  y  as 

“U  *  /ux»  Y)  x 

parameters.  These  graphs  are  contained  in  Schubel ,  Carter  et  al., 

(1978)  and  some  are  included  in  Chapter  V  of  this  report  as  examples. 

The  lateral  dimensions  of  the  plume  are  determined  by  taking  the 

second  moment,  y»  of  the  concentration  distribution  of  Equation  2.4. 

The  second  moment  can  also  be  described  as  a  function  of  x*.  and 

x 

y  and  has  the  value: 

F  (x>t)  =  ^  F  (x*.  “Vux,  y)  (2.7) 

where 

/  t*  exp  -  (rf)  f—~  t  )  exp  -  (yt*)  dt* 

F(x*.  U/Ux>  y)  =  ^ - 7— - - -  (2.8) 

f  r  exp  -  exP  "  ^t*)  dt* 

Again,  Equation  2.7  is  normalized  with  respect  to  the  second  moment  at 
the  plume  front  to  obtain: 


The  lateral  dimensions  of  the  plume  are  determined  from  another  set  of 
graphs  in  Schubel,  Carter  et  al.,  (1978).  A  few  examples  are  shown 
in  the  example  calculation  in  Chapter  V  of  this  report. 


C(x,0,t) 

C(uxt,o,t 
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This  model  is  particulary  applicable  to  dredge  disposal  in  a 

shallow,  wide  estuary.  The  assumption  of  K  =  K  is  only  valid  in  an 

x  y 

area  where  there  is  not  a  strong  primary  flow  velocity,  u  .  The  ver- 
tically  averaged  suspended  solids  concentration  is  suitable  for  a  shal¬ 
low  disposal  area.  The  assumption  of  a  vertical  line  source  is  also 
typical  of  the  normal  mode  of  dredge  disposal  in  an  estuarine  environ¬ 
ment  (Barnard,  1978). 


Wechsler  and  Cogley  Model 

The  model  developed  by  Wechsler  and  Cogley  (1977)  is  for  predic¬ 
tion  of  downstream  concentration  of  suspended  sediment  in  waters 
characterized  by  unidirectional,  steady  flow,  infinite  width,  constant 
depth  and  infinite  length.  The  initial  differential  equation  for 
describing  the  suspended  solids  concentration  at  any  point  downstream 
of  the  dredge  discharge  is: 

h  (uxc)  +  h  (fWf<“>dW)  -  h  (Kx  §)  -  57  (K,  If) 

-  h  (Kz  If)' 0  <2-10> 

where  x,y  and  z  represent  the  longitudinal,  lateral  and  vertical  co¬ 
ordinates,  respectively;  ux  is  the  mean  current  velocity  in  the  x 
direction;  C  is  the  suspended  sediment  concentration;  W  is  the  settl¬ 
ing  velocity;  K  ,  K  ,  and  K  are  the  eddy  diffusion  coefficients  in 
X  y  z 

the  x,  y,  and  z  directions;'  and  f(W)  is  the  settling  velocity  frequency 
distribution.  The  first  term  in  Equation  2.10  describes  downstream 
advection,  the  second  term  describes  vertical  sedimentation, whi le  the 
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last  three  terms  describe  eddy  diffusion  in  the  x,  y  and  z  directions, 
respectively. 

Several  simplifying  assumptions  are  made  to  make  the  model  useful: 
1)  the  eddy  diffusion  in  the  downstream  direction  is  negligible  com¬ 
pared  to  the  other  diffusion  and  transport  terms,  therefore. 


3  (  3C  \ 

3x  Kx  dx)  =  ^  the  e<ic*-y  dlf^us10n  the  vertical  direction  can  be 

related  to  the  vertical  position  in  the  flow  by,  Kz  =  0.02  uxz  ^1-jA  ; 
3)  the  eddy  diffusion  in  the  lateral  direction  is  given  by 
K  -  2.2  (K  1  :  and  4)  for  non-flocculant  sediment,  the  settling  term 

y  Z  mu  A 

can  be  described  by  W  3C/3z,  and  solving  the  model  for  each  sediment 
size  fraction  and  superimposing  the  results  for  the  final  solution. 

The  resulting  equation  is: 


u  5£  +  u?c  _a_ 
ux  3x-  3z  3Z 


(0.02  uxz  0-§)  |f)  -  fj-  (2.2  (K2)max  |f)  = 


0 

(2.11) 


Equation  2.11  is  solved  using  the  finite  difference  method  for  the 
downstream  and  vertical  directions  and  an  analytical  solution  involv¬ 
ing  the  "error  function"  for  the  lateral  direction.  It  is  assumed  the 
source  is  a  vertical  line  source,  continuously  emitting  sediment  at  a 
given  strength  per  unit  height.  This  source  strength  is  converted  to 
a  concentration  by  assuming  the  sediment  is  initially  concentrated  in 
a  vertical  column  of  width,  b,  which  is  small  relative  to  the  depth,  D. 
The  upstream  boundary  condition  is  then.C  =  CQ  at  x  -  0,  |y|  <  b, 

z  <_  D.  The  surface  boundary  condition  specifies  no  net  flux  of  material 

3C 

across  the  surface,  or  gy  +  WC  =  T  ,e  bottom  boundary  condition 


assumes  all  material  settling  to  the  bottom  remains,  with  no  re-entrain¬ 
ment,  or  Kz  =  0* 

The  model  solution  is  contained  in  a  computer  program  which  is 
described  in  Wechsler  and  Cogley  (1977).  The  inputs  to  the  program 
are  mean  current  velocity,  mean  stream  depth,  settling  velocity  dis¬ 
tribution  (given  as  any  number  of  sediment  fractions  and  their  corre¬ 
sponding  concentration  and  settling  velocity)  and  three  computational 
parameters.  The  output  consists  of  1)  a  section  showing  the  vertical 
distribution  of  sediment  downstream  for  each  sediment  fraction  (without 
lateral  spreading);  2)  the  summation  of  the  vertical  distributions  for 
all  size  fractions;  3)  the  lateral  spreading  coefficients;  and  4)  hori¬ 
zontal  slices  through  the  three-dimensional  plume  at  five  pre-selected 
depths  showing  the  concentration  distribution  at  each  depth. 

The  assumption  of  a  vertical  line  source  of  width,  b,  which  is 

less  than  the  total  depth,  D,  is  applicable  to  open  water  discharge 

of  dredged  material.  It  is  less  applicable  to  a  plume  resulting  from 

land  runoff  since  the  plume  source  tends  to  be  wide  with  respect  to 

the  depth.  The  assumptions  concerning  the  eddy  diffusivities,  K  ,  K 

X  y 

and  Kz>  are  suitable  for  describing  a  plume  developing  in  a  river  or 
an  estuary  with  a  strong  current  flow. 

Convection  -  Dispersion  Equation 

The  basic  equation  describing  convection  and  dispersion  of  dis¬ 
solved  matter  or  suspended  particles  is  based  on  the  principle  of  con¬ 
servation  of  mass.  For  a  conservative  substance,  the  principle  of  con¬ 
servation  of  mass  can  be  stated  (Sayre,  1968): 


.1 

>4 

[i 
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Rate  of  change" 

Rate  of  change  of 

Rate  of  change  of 

of  mass  in 

= 

mass  in  control 

+ 

mass  in  control 

_control  volume 

volume  due  to 

volume  due  to 

convection 

diffusi on 

9C 

at 

= 

-  u.  |£ 
l  9xi 

+ 

where  e.  is  the  diffusion  coefficient  in  the  i'th  direction  and  all 
other  terms  are  described  previously.  For  laminar  flow,  =  e^,  the 

coefficient  of  molecular  diffusion.  For  turbulent  flow,  e.  =  eT  +  eM, 
where  e-j-  is  the  coefficient  of  turbulent  diffusion.  In  Fickian  diffu¬ 
sion  theory,  it  is  assumed  that  dispersion  resulting  from  turbulent 
open-channel  flow  is  exactly  analogous  to  dispersion  from  molecular 
diffusion.  The  dispersion  coefficients  in  the  x,  y,  and  z  directions 

are  assumed  to  be  constants,  given  by  K  ,  K  and  K  .  The  resulting 

x  y  z 

equation,  expressed  in  Cartesian  coordinates  is: 


9C  .  „  9C  ,  „  9C  ,  9C  _  „  92C  .  „  92C  ,  „  92C 
3t  Ux  97  +  uy  97  +  Uz  Ti  -  Kx  aF  +  Ky  9^  +  Kz  9?? 


(2.13) 


The  solution  of  Equation  2-.  13  depends  on  the  values  of  K  ,  K  and 

x  y 

«z.  Various  authors  have  arr.ved  at  equations  to  approximate  the 
values  of  the  dispersion  coefficients  (K)  in  the  longitudinal  (x), 
lateral  (y),  and  vertical  (z)  directions. 


Longitudinal  Dispersion  Coefficient 
The  first  discussion  of  dispersion  in  turbulent  flow  was  by  Taylor 
(1954)  for  dispersion  in  a  long,  straight,  circular  pipe.  Taylor  found 
the  dispersion  coefficient  to  be: 
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Kx  =  10.06  a  U*  (2.14) 

where  a  is  the  pipe  radius  and  U*  is  the  shear  velocity.  The  shear 
velocity  can  be  calculated  by  U*  = Vx  /p,  where  is  the  shear  stress 
at  the  wall  of  the  pipe  and  p  is  the  fluid  density. 

Elder  (1959)  obtained  an  expression  for  Kx  in  two-dimensional  open- 
channel  flow: 

Kx  =  a  D  U*  .(2.15) 

where  a  =  5.93,  D  is  the  mean  depth  and  U*  is,  again,  the  shear  velocity, 

calculated  as  U*  =  Vtq/p  =  Vg  D  Se  where  tq  is  the  shear  stress 
at  the  bottom,  g  is  the  acceleration  of  gravity  and  Sg  is  the  energy 
slope.  Elder's  expression  is  for  infinitely  wide  channels,  meaning  no 
lateral  velocity  or  concentration  gradients,  and  a  logarithmic  ver¬ 
tical  velocity  distribution.  Longitudinal  dispersion,  therefore,  is 
a  result  of  differential  convection  in  the  vertical  direction  and 
turbulent  diffusion. 

Yotsukura  and  Fiering  (1964)  applied  Taylor's  solution  method  to 
open  channels  and  used  a  computer  solution  to  obtain  values  of  a  vary¬ 
ing  from  9  to  13  as  the  ratio  of  ux/U*  varied  from  14.5,  indicating  a 
rough  channel  boundary,  to  20,  indicating  a  smooth  channel  boundary. 

Thackston  and  Krenkel  (1967)  included  the  term  ux/U*  in  the  dis¬ 
persion  equation,  resulting  in: 

‘1/4 

Kx  =  ctDU*Vw  (2J6) 

where  a  has  the  value  5.82  or  7.25.  The  value  ux/U*  is  a  dimensionless 
measure  of  the  bottom  roughness;  larger  values  meaning  smoother  bottoms. 
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Thackston  and  Krenkel  are  careful  to  point  out,  however,  that  Equation 
2.16,  as  well  as  all  of  the  previously  mentioned  equations,  does  not 
apply  in  areas  where  there  is  appreciable  lateral  velocity  variation. 

In  such  a  case,  the  authors  state  that  K  will  be  much  larger  than 

A 

calculated  by  Equation  2.16,  and  recommend  in  situ  measurement  of 
Kx>  Since  natural  streams  have  a  significant  lateral  velocity  profile, 
none  of  the  preceeding  equations  and  a  coefficients  are  directly 
applicable. 

Fischer  (1966)  showed  that  the  dispersion  of  a  slug  of  material 
injected  into  a  natural  stream  is  divided  into  two  distinct  phases; 

1)  the  convective  period,  in  which  the  material  diffuses  laterally  and 
longitudinally  until  the  material  is  completely  distributed  across  the 
channel,  and  2)  the  diffusive  period  (called  the  Taylor  period),  in 
which  the  lateral  concentration  gradient  is  small.  The  convective 
period  is  characterized  by  a  highly  skewed  longitudinal  concentration 
profile;  the  downstream  face  being  blunt  and  the  upstream  tail  being 
long.  The  above  equations  for  Kx  are  not  applicable  to  the  convective 
period.  The  Taylor  period  is  characterized  by  a  more  nearly  Gaussian 
longitudinal  concentration  profile.  The  above  equations  are  applicable 
with  the  restrictions  mentioned,  to  the  Taylor  period.  The  criterion 

for  determining  if  dispersion  of  a  material  is  in  the  convective  period 

* 

or  the  Taylor  period  is'  (Fischer,  1966): 
l2  Uy 

L  >  1.8-  ~  (2.17) 

where  L  is  the  distance  downstream  from  the  source  of  the  matnrial;  Z 
is  the  characteristic  cross-sectional  length,  described  as  the  distance 


14 


dispersion  coefficient,  taken  as  Ky  =  0.23  DU*  by  Fischer  (1967). 
Equation  2.18  can  be  solved  for  any  stream  after  measuring  the  energy 
slope,  Sg,  the  cross-sectional  geometry  and  the  cross-sectional  vel¬ 
ocity  distribution  of  a  "typical"  cross-section.  Fischer  (1967)  solved 
Equation  2.18  with  the  use  of  a  computer  for  several  laboratory  flumes 
and  related  the  resulting  Kx  values  back  to  Equation  2.15  and  found 
values  of  a  ranging  from  5  to  16.  The  higher  values  of  a  were  for 
flumes  with  sloping  sides  rather  than  perpendicular  sides.  Again,  the 
lateral  velocity  currents  set  up  by  the  sloping  sides  of  natural  streams 
give  problems  in  predicting  a,  so  the  more  simple  Equation  2.15  can  not 
be  used. 


Liu  (1977)  used  Equation  2.18,  since  it  correctly  describes  the 
prime  mechanism  of  dispersion  in  natural  streams,  to  develop  an  ex¬ 
pression  for  Kx  which  is  much  easier  to  calculate: 


u  2  B3 

Kx ' p  "for 


(2.20) 


where  (Liu  1978  ) , 

B  =  O.502 

and  QB  is  the  river  discharge.  The  new  coefficient,  8,  is  an  easier 
coefficient  to  use  than  a,  since  8  does  not  depend  on  stream  morpho¬ 
metry  but  on  the  dimensionless  bottom  roughness,  a  value  more  easily 
estimated,  based  on  existing  data  for  Kx  in  streams  and  the  value  of 
K  predicted  by  Equation  2.20,  K  can  be  predicted  to  within  a  factor 

A  A 

of  six  by  Equation  2.20.  This  is  better  than  any  other  of  the  simple 
methods  described  for  predicting  the  longitudinal  dispersion  coeffi¬ 


cient. 
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Lateral  Dispersion  Coefficient 

Elder  (1959)  proposed  the  equation  for  predicting  the  lateral  dis¬ 
persion  coefficient,  Ky: 

Ky  =  <P  D  U*  (2.22) 

where  <j>  is  equal  to  0.23.  The  value  of  <J>  =  0.23  was  obtained  by  ex¬ 
periment  in  long,  wide  laboratory  flumes. 

Many  authors  have  since  investigated  the  value  of  <J>  in  both  labor¬ 
atory  flumes  and  natural  streams.  Sayre  and  Chang  (1968)  reported 
<J>  =  0.17  in  a  straight  laboratory  flume.  Yotsukura  and  Cobb  (1972) 
report  values  of  <f>  for  natural  streams  and  irrigation  canals  varying 
from  0.22  to  0.65,  with  most  values  being  near  0.3.  Other  reported 
values  of  <j>  range  from  0.17  to  0.72.  The  higher  values  for  c>  are  all 
for  very  fast  rivers.  The  conclusions  drawn  are  that;  1)  the  form  of 
Equation  2.22  is  correct  for  predicting  K^,  but  <|>  may  vary,  and  2)  ap¬ 
plication  of  Fickian  theory  to  lateral  dispersion  is  correct  as  long  as 
there  are  no  appreciable  lateral  currents  in  the  stream. 

Okoye  (1970)  refined  the  determination  of  <J>  somewhat  by  use  of  the 
aspect  ratio,  X  =  D/B,  the  ratio  of  the  stream  depth  to  stream  width. 

It  was  found  that  <f>  decreased  from  0.24  to  0.093  as  X  increased  from 
0.015  to  0.200. 

The  effect  of  bends  in  the  channel  on  is  significant.  Yotsukura 
and  Sayre  (1976)  reported  that  $  varies  from  0.1  to  0.2  for  straight 
channels,  ranging  in  size  from  laboratory  flumes  to  medium  size  irri¬ 
gation  chanals;  4>  varies  from  0.6  to  10  in  the  Missouri  River,  and  <t> 
varies  from  0.5  to  2.5  in  curved  laboratory  flumes.  Fischer  (1968) 
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from  the  point  of  maximum  surface  velocity  to  the  far  bank;  r  is  the 
hydraulic  radius  and  ux/U*  is  as  previously  defined.  If  L  is  greater 
than  the  right  hand  side  of  Equation  2.17,  then  the  Taylor  period  has 
been  reached. 

Working  with  natural  streams,  Fischer  (1967)  found  that  longitud¬ 
inal  dispersion  was  a  result  of  the  combination  of  two  effects; 

1)  variable  lateral  convective  velocities  and  2)  concentration  grad¬ 
ients  giving  rise  to  lateral  diffusion  of  material.  The  effect  of  the 
lateral  diffusion  is  to  dampen  the  dispersion  caused  by  the  different¬ 
ial  lateral  convective  velocities.  This  mechanism  for  dispersion  in 
natural  streams  is  in  contrast  to  the  mechanism  proposed  by  Elder  (1954) 
and  used  by  the  other  authors,  in  that  dispersion  is  caused  by  lateral 
velocity  gradients  as  opposed  to  vertical  velocity  gradients. 

Using  -this  mechanism,  Fischer  (1967)  found  an  equation  for  the 
longitudinal  dispersion  coefficient  in  the  Taylor  period: 


*  / 


q’(y)  dy 


/  V' 


q'(y)  dy  (2.18) 


where 


q'(y)  = 


u'(z,y)  dz 


(2-19) 


and  q'(y)  is  described  as  the  discharge  per  unit  width;  u’  is  the  de¬ 
viation  of  the  local  mean  velocity,  U,  from  the  cross-sectional  mean 

velocity,  u  ,(U  -  uv  -  u');  B  is  the  stream  width;  and  K  is  the  lateral 
xx  y 
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Lateral  Dispersion  Coefficient 

Elder  (1959)  proposed  the  equation  for  predicting  the  lateral  dis¬ 
persion  coefficient,  Ky: 

Ky  -  0>  D  U*  (2.22) 

where  <j>  is  equal  to  0.23.  The  value  of  =  0.23  was  obtained  by  ex¬ 
periment  in  long,  wide  laboratory  flumes. 

Many  authors  have  since  investigated  the  value  of  <j>  in  both  labor¬ 
atory  flumes  and  natural  streams.  Sayre  and  Chang  (1968)  reported 
<f>  =  0.17  in  a  straight  laboratory  flume.  Yotsukura  and  Cobb  (1972) 
report  values  of  <j>  for  natural  streams  and  irrigation  canals  varying 
from  0.22  to  0.65,  with  most  values  being  near  0.3.  Other  reported 
values  of  4>  range  from  0.17  to  0.72.  The  higher  values  for  <J>  are  all 
for  very  fast  rivers.  The  conclusions  drawn  are  that;  1)  the  form  of 
Equation  2.22  is  correct  for  predicting  Ky,  but  <J>  may  vary,  and  2)  ap¬ 
plication  of  Fickian  theory  to  lateral  dispersion  is  correct  as  long  as 
there  are  no  appreciable  lateral  currents  in  the  stream. 

Okoye  (1970)  refined  the  determination  of  <j>  somewhat  by  use  of  the 
aspect  ratio,  X  =  D/B,  the  ratio  of  the  stream  depth  to  stream  width. 

It  was  found  that  <J>  decreased  from  0.24  to  0.093  as  X  increased  from 
0.015  to  0.200. 

The  effect  of  bends  in  the  channel  on  Ky  is  significant.  Yotsukura 
and  Sayre  (1976)  reported  that  4>  varies  from  0.1  to  0.2  for  straight 
channels,  ranging  in  size  from  laboratory  flumes  to  medium  size  irri¬ 
gation  chanals;  <f>  varies  from  0.6  to  10  in  the  Missouri  River,  and  4> 
varies  from  0.5  to  2.5  in  curved  laboratory  flumes.  Fischer  (1968) 


reports  that  higher  values  of  4>  are  also  found  near  the  banks  of 
rivers. 


Vertical  Dispersion  Coefficient 

Very  little  experimental  work  has  been  done  on  the  vertical  dis¬ 
persion  coefficient,  Kz.  Jobson  and  Sayre  (1970)  reported  a  value 
for  marked  fluid  particles  of: 

Kz  =  «U*z  (l  -  §)  (2.23) 

for  a  logarithmic  vertical  velocity  distribution,  k  is  the  von  Karman 
coefficient,  which  is  shown,  experimentally,  to  be  approximately  =  0.4 
(Tennekes  and  Lumley  1972  ).  Equation  2.23  agrees  with  experimental 
data  fairly  closely. 


Water  Quality  Criteria 

The  federal  water  quality  criterion  for  turbidity  and  suspended 
solids  is  based  on  protection  of  freshwater  fish  and  other  aquatic 
life  (Water  Quality  Criteria  1976  ).  The  criterion  is  stated: 
"settleable  and  suspended  solids  should  not  reduce  the  depth  of  the 
compensation  point  for  photosynthetic  activity  by  more  than  10  percent 
from  the  seasonally  established  norm  for  aquatic  life." 

Turbidity  and  suspended  solids  have  several  effects  on  fish  and 
other  aquatic  organisms.  Deposited  sediments  can  damage  invertebrate 
populations  and  cover  gravel  spawning  areas.  Silt  attached  to  eggs 
may  inhibit  oxygen  transfer  and  so  increase  mortality.  Suspended  sol¬ 
ids  may  act  directly  on  fish  by  either  killing  them  or  inhibiting 
their  growth,  and  by  reducing  the  availability  of  food.  Suspended 
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solids  reduce  light  penetration  which  causes  a  reduction  in  the  depth 
of  the  photic  zone.  This  reduced  photic  zone  may  lead  to  a  reduction 
in  primary  production  which  leads  to  a  decrease  in  the  amount  of  food 
for  fish.  Turbidity  also  interferes  with  aesthetic  enjoyment  of  water¬ 
ways. 

The  Iowa  Water  Quality  Standard  (1977)  for  surface  water  states; 
"the  turbidity  of  the  receiving  water  shall  not  be  increased  by  more 
than  25  Nephelometric  turbidity  units  by  any  point  source  discharge." 

The  criterion  shall  apply  after  an  appropriate  mixing  zone.  The  mix¬ 
ing  zone  is  the  area  of  diffusion  of  an  effluent  in  the  receiving  water. 
In  all  cases,  the  mixing  zone  should  be  as  small  as  practicable  and 
not  include  more  than  25  percent  of  the  cross-sectional  area. 
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CHAPTER  III 
MODEL  DEVELOPMENT 


The  distribution  of  sediment  in  the  water  column  is  governed  by 
the  equation: 


9C  +  u.  JC_ 
at  i  axi 


3xi 


V  3C  ' 
Ki  ax. 


(3.1) 


where  ui  refers  to  the  fluid  velocity  in  a  rectangular  coordinate 
system  (xi ),  C  refers  to  the  concentration  (mass  per  unit  volume)  of 
sediment  suspended  in  the  water  column,  and  refers  to  the  dispersion 
coefficient  in  the  i'th  direction  (Sayre  1968).  Equation  (3.1)  can  be 
rewritten  as: 

/  \  /  \  /  \ 

ac 


sc  ,  u  ac  ,  u  ac  ,  u.  ac  .  a  (  v 
at  ax  y  ay  z  az  ax  x  ax 


+  9_(K  W  a 
sy  \Ky  ay  / 


K  I?) (3.2) 


3z  y'z  3z/ 

where  x,  y  and  z  refer  to  the  longitudinal,  lateral  and  vertical 
directions,  respectively. 

The  solids  in  the  plume  are  not  uniform,  but  consist  of  various 
size  particles,  each  with  a  distinct  settling  velocity.  The  vertical 
velocity  of  a  particle  (uz)  can  be  divided  into  two  fractions,  its 
natural  settling  velocity  in  quiescent  water,  w^ ,  and  the  velocity  of 
the  water  in  the  z  direction,  w.  Incorporating  these  into  Equation 
(3.2)  gives: 
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ft  Zci  +  ux  h  Zci  +  uy  |r  Eci  +  K-  +  w>  If  Zci  = 


tx  (K*ljEci)  +  ly  (Ky  fy  Eci)  *h  (K2flE<:i)  <3-3> 

These  plumes  develop  along  a  shore  of  the  river  where  the  water  is 

shallow;  therefore  vertically  averaged  solids  concentrations  will  be 

calculated.  The  necessary  assumptions  are  that  u  ,  u  ,  K  and  K  are 

x  y  x  y 

depth  independent,  w  =  0  and  there  can  be  no  flux  of  suspended  material 

HP 

across  the  surface  of  the  river  (-  Kz  +  uzC=  0  at  z  =  0).  With 
these  assumptions,  Equation  (3.3)  can  be  integrated  to  obtain  (Schubel , 
et  al.,  1978): 


3C  .  „  3C  .  8C_3  L  3CV3_/£  3C\  WC 
8t  x  3x  Uy  3y  3x  lx  3x/  3y  lKy  3y/'  D 


(3.5) 

(3.6) 

(3.7) 


C  is  defined  as  the  mean  suspended  solids  concentration  and  W  is  the 
mean  settling  velocity  Of  the  particles.  D  is  the  average  depth  of 
the  water  containing  the  plume. 

In  a  river,  the  following  additional  assumptions  can  be  made  to 
further  simplify  Equation  (3.4)  (Sayre,  1973). 


20 


u  =  constant 

A 

D  =  depth  =  constant 
“y-° 


Ky  =  constant 

and  K  «  K  so 
x  y 


i _L  3C ' \  _ 
3x  l  x  3x J 


0 


(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 


Incorporating  these  assumptions  into  Equation  (3.4)  and  assuming  steady 
state  (3C/3t  =  0)  gives: 


3C  =  Kw  3_C 
3x  "  uj  3yr 


WC 

uxD 


(3.13) 


The  solution  to  this  equation  can  be  written  (Sayre,  1979)  as: 

Wx 


uxD 


C  (y,x)  =  C'(y,x)  exp 
where  C’(y,x)  is  the  solution  to  the  diffusion  equation: 


(3.14) 


3C._Kw.3_C 
ux  dy2 


(3.15) 


For  the  case  of  a  continuous  point  source  of  flow,QQ  and  solids 
concentration, CQ,  the  solution  to  Equation  (3.15)  is  (Sayre,  1973): 


C‘ (y,x)  = 


Q0C 


o^o 


1 


UxD  2\|  7T  Ky  x/u, 


exp 


2 

y  u 
J  x 

4  K  x  , 

•-  y  -1 


(3.16) 


This  equation  has  the  form  of  a  normal  probability  function  with  var¬ 
iance,  a  2  =  2K  x/u  .  .Substituting  this  into  Equation  (3.16)  gives: 

y  y  x 


C  (y.x)  = 


Vo 


exp 


[^*1 


(3.17) 


X 


Equation  (3.17)  describes  the  plume  resulting  from  a  continuous 
point  source.  Water  running  off  from  a  shore  can  better  be  described 
by  a  line  source  of  width, b,  perpendicular  to  the  shoreline.  Equation 
(3.17)  can  be  modified  to  describe  a  line  source  by  the  method  used  by 
Sayre  (1973,  1979).  The  resulting  expression  is: 

B  r 

C'  (y.x)  =  /  C'(y\0)  —J -  exp  -  -4^-  dy'  (3.18) 

o  oy  L  y  _ 

where  B  is  the  width  of  the  river  and  y’  is  a  dummy  variable  describing 
any  point  within  the  source  width.  The  initial  conditions  for  this 
line  source  are: 


C'  (y'  ,0)  =  * 

C’  (y',0)  =  0 


0  <  y'  <  b 


b  <  y'  <  B 


(3.19) 


where  Qb  =  ‘i^Db  and  is  the  portion  of  the  river  flow  passing  through 
the  source  width, b.  Incorporating  these  into  Equation  (3.18)  and  sub¬ 
stituting  the  standard  normalized  variable  s  =  gives 

r  y/oy  °y  1 

nr  r  •>  i—  i 


C'  (y.x)  = 


0  0  1 


b  vfz^T 


si 

2  J 


(3.20) 


which  is  in  the  form  of  the  cumulative  normal  distribution  function. 


The  solution  to  Equation  (3.20)  is 


F  (J-)  .  F  (tb) 

L  y  yJ 


(3.21) 


C '  (y.x)  =  F  (J-)  -  F  (3.21) 

yb  L  y  yJ 

where  the  value  o*  F  (*)  can  be  obtained  from  a  cumulative  normal  dis¬ 


tribution  table,  such  as  the  one  included  as  Appendix  C. 


The  suspended  solids  plume  described  by  Equation  (3.21)  includes 
no  effects  from  the  side  banks  of  the  river.  It  is  assumed  that  the 
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channel  banks  act  as  reflecting  barriers.  Including  the  effects  of 
reflection  from  the  near  side  bank,  the  equation  becomes: 

C*  (y.x)  =  ^  F  (*-^)  -  F  (3.22) 

%  L  y  y  J 

This  equation  is  not  applicable  if  the  suspended  solids  plume  disperses 
in  the  lateral  direction  enough  to  reflect  from  the  far  shoreline.  An 
exact  solution  is  presented  by  Sayre  (1969). 

Substituting  Equation  (3.22)  back  into  Equation  (3.14)  gives  the 


final  solution. 


Qoco 


c  (y,x)  =  F  ' F  {X^r)  exp  ■  S  (3-23) 

is  the  initial  suspended  solids  concentration  at  the  source.  By 


dividing  both  sides  of  Equation  (3.23)  by  the  initial  concentration,  the 


model  can  be  written: 


C  (y.x)  ■  F  (^>  -  F  exp  -  jji  (3.24) 

and  the  right  hand  side  can  be  solved  independent  of  the  source  concen¬ 
tration. 

The  parameters  that  are  necessary  to  solve  the  model  are  the 
source  width,  b;  the  mean  depth,  D;  the  mean  downstream  velocity,  ux; 
the  lateral  dispersion  coefficient,  and  the  terminal  settling 
velocity  of  the  suspended  particle,  W.  Values  of  downstream  distance, 


i 


x,  are  chosen  and  the  lateral  extent  of  the  plume  is  calculated  by  vary¬ 
ing  the  value  of  y/o  ,  and  hence,  y.  An  additional  advantage  is  that 
the  model  can  be  solved  several  times  for  size  fractions  with  different 
terminal  settling  velocities  and  the  several  solutions  summed  for  the 
final  solution, due  to  the  principle  of  superposition  for  linear  dif¬ 
ferential  equations. 

The  model  can  be  programmed  for  solution  with  a  programmable  cal¬ 
culator.  One  program  for  a  Hewlett-Packard  29C  is  included  as 
Appendix  A. 
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CHAPTER  IV 

FIELD  PROCEDURES  AND  RESULTS 


Suspended  Solids/Turbidity  Relationships 
It  was  intended  to  use  a  continuous  flow  turbidity  monitoring 
device  to  sample  the  plume.  Discrete  samples  were  also  to  be  taken 
and  analyzed  for  suspended  solids  concentration.  With  this  data,  a 
correlation  could  be  developed  to  translate  the  continuous  flow  tur¬ 
bidity  data  into  suspended  solids,  which  was  necessary  for  input  into 
the  model.  To  this  end,  experiments  were  carried  out  in  the  laboratory 
to  develop  correlations  for  three  distinct  types  of  particles,  sand, 
laboratory  grade  colloidal  kaolin  clay*  and  Iowa  River  mud,  a  mixture 
of  silt  and  clay. 

Turbidity  was  measured  nephelometrically  with  a  Turner  Model  111 
Fluoiometer  equipped  with  a  flow-through  door.  A  2A  secondary  filter 
was  used  with  no  primary  filter.  The  sample  of  turbid  water  was  con¬ 
tained  in  a  1000  ml  Erlenmeyer  flask  and  was  continuously  mixed  with  a 
magnetic  stirrer  and  stir  bar.  The  sample  was  withdrawn  from  the  flask, 
drawn  through  the  fluorometer  at  approximately  1.2  1/min,  and  returned 
to  the  flask.  When  a  steady  turbidity  reading  was  obtained,  a  sample 
was  collected  from  the  pump  discharge  and  analyzed  for  suspended  solids. 
Flow  was  downflow  through  the  fluorometer.  The  material  in  the  flask 

*  Fisher  Scientific  Co.,  Fair  Lawn,  N.J. ,  Laboratory  Grade  Colloidal 
Kaolin  Powder. 


was  then  diluted  and  the  procedure  repeated.  Figure  4-1  shows  the 
equipment  used. 
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The  procedure  was  repeated  for  each  of  the  three  types  of  mater¬ 
ials,  sand,  kaolin,  and  river  mud.  The  relationship  between  suspended 
solids  and  turbidity  for  each  of  these  materials  is  shown  in  Figure 
4-2.  It  can  be  seen  that  although  there  are  great  differences  in  the 
suspended  solids  concentration  necessary  to  produce  a  certain  turbidity, 
each  material  exhibits  a  distinct  relationship  between  suspended  solids 
and  turbidity.  The  clay  particles  are  smaller  and  more  numerous  per 
unit  mass  and  therefore  scatter  light  to  a  greater  degree  than  Iowa 
River  mud  or  sand. 

It  was  felt  the  suspended  solids  in  the  plume  resulting  from 
disposal  operations  would  exhibit  this  same  phenomenon.  It  was  there¬ 
fore  decided  to  measure  turbidity  continuously  in  transects  across 
the  plume  and  take  enough  discrete  samples  for  suspended  solids  analysis 
to  describe  the  relationship  between  the  two  parameters. 

Field  Sampling 

Three  of  the  four  dredging  operations  by  the  U.5.  Army  Corps  of 
Engineers,  Rock  Island  District,  on  the  Mississippi  River  in  1978  were 
monitored.  Dredging  operations  monitored  were  near  Hannibal,  Missouri, 
river  mile  313.5,  on  October  16  and  17;  near  Keithsburg,  Illinois, 
river  mile  425.8,  on  October  25;  and  at  Rock  Island,  Illinois,  river 
mile  482.0  on  October  28.  These  three  sites  are  shown  in  Figure  <:i-3. 

At  the  Hannibal  site,  18,800  cubic  yards  of  sediment  were  dredged.  At 
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Figure  4-2.  Turbidity  vs.  suspended  solids  relationships  for  sand 
colloidal  kaolin  and  Iowa  River  mud. 
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Keithsburg,  11,166  cubic  yards  of  sediment  were  dredged  in  16.0  hours 
and  at  Rock  Island,  11,596  cubic  yards  of  sediment  were  dredged  in 
18.58  hours.  Each  dredging  disposal  operation  was  unique  with  respect 
to  the  resulting  turbidity  plume  generated. 

The  dredge  spoil  at  the  Hannibal  site  was  discharged  to  nearby 
Armstrong  Island,  shown  in  Figure  4-4.  This  island  is  approximately 
1.7  miles  long  and  0.3  miles  wide  at  its  widest  point.  It  also  has  a 
large  inland  depression  and  lake.  There  was  no  runoff  from  this 
island  during  the  dredging.  Much  of  the  discharged  water  was  assumed 
to  be  percolating, with  the  rest  ponding  in  depressions  on  the  island. 
Samples  of  the  discharged  water  and  of  the  ponded  water  were  collected 
for  size  analysis  of  the  suspended  solids  for  the  purposes  of  comparison. 

The  dredge  spoil  at  the  Keithsburg  site  was  discharged  to  Willow 
Bar  Is  1  and ,_  adjacent  to  the  dredge  cut,  see  Figure  4-5.  Willow  Bar 
Island  is  approximately  2500  feet  long  and  400  feet  wide  and  gently 
slopes  away  from  the  main  channel  of  the  river.  Consequently,  there 
was  a  return  water  flow  to  the  back  side  of  the  island.  Several  points 
of  entry  were  noted  but  only  the  area  downstream  from  the  major  runoff 
point  was  monitored. 

The  possible  lateral  and  longitudinal  dimensions  of  the  turbidity 
plume  were  estimated  from  surface  debris  washed  into  the  river  with 
the  runoff  flow.  A  system  of  shore  markers  and  in-stream  buoys  was 
laid  out  to  act  as  location  markers  so  that  the  dimensions  of  the 
plume  could  be  accurately  determined.  A  geodimeter  (distance  meter), 
Hewlett-Packard  Model  3800B,  was  used  to  measure  the  distance  of  each 
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of  the  shore  markers  and  buoys  from  the  source  of  the  plume.  Each  buoy 
was  placed  so  that  it  was  roughly  perpendicular  from  the  shore  line. 

With  this  information,  it  was  possible  to  calculate  all  of  the  distances 
in  the  grid  system  of  markers  and  buoys.  The  grid  system  was  approx¬ 
imately  180  m.  long  by  80  m.  wide. 

Sampling  of  the  plume  was  done  by  traversing  the  grid  system  in  a 
serpentine  fashion  in  a  slow  moving  boat.  Samples  were  drawn  through 
the  fluorometer  continuously  with  a  Masterflex  Model  7545,  Variable 
Speed  Drive  pump.  The  pump  was  equipped  with  a  number  7017  head  and 
used  0.225  in.  I.D.  by  0.3900  in.O.D.  Tygon tubing.  Samples  were  drawn 
at  a  rate  of  approximately  0.6  1/min  and  had  an  approximate  residence 
time  of  0.4  minutes  in  the  tubing.  The  boat  was  estimated  to  be  moving 
at  l  m/sec  so  the  boat  had  moved  approximately  20  m.  between  the  time 
the  sample -was  removed  from  the  water  column  and  the  time  the  turbidity 
was  read  and  the  sample  collected  for  suspended  solids  analysis.  A 
YSI  Model  81A  recorder  with  a  30  in/hr  chart  speed  gear  was  attached 
to  the  fluorometer  to  continuously  record  the  turbidities.  100  ml 
discrete  samples  were  taken  from  the  pump  discharge  for  calibration  of 
the  turbidity  vs.  suspended  solids  relationship. 

The  plume  was  sampled  at  three  depths;  top,  middle  and  bottom. 
During  the  sampling,  the  fluorometer  became  inoperable.  It  began 
showing  relatively  constant  turbidity  readings  at  all  points  in  the 
plume.  It  was  also  giving  an  abnormally  high  reading  for  the  turbidity, 
around  500  to  900  NTU.  Normal  turbidity  readings  were  all  less  than 
100  NTU.  Consequently,  the  continuous  output  was  not  used. 
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A  i<  itional  water  samples  were  collected  for  size  analysis  of  the 
susp.---.ded  sediment  at  the  head  of  the  plume,  and  at  the  discharge  point 
into  the  river.  Size  analysis  included  visual  accumulation  tube  (VA 
tube)  as  well  as  micropipette  measurements  for  coarse  and  fine  graded 
materials.  The  velocity  of  the  water  flowing  into  the  river  was  suf¬ 
ficient  to  erode  the  shoreline  of  the  island.  A  channel  was  cut  into 
the  shoreline  approximately  three  feet  wide  at  the  mouth  and  extending 
approximately  fifteen  feet  inland.  Since  this  material  was  forming 
the  plume,  a  sample  of  this  soil  was  collected.  At  the  point  where 
the  flow  entered  the  river,  a  sand  bar  was  built  up  during  the  course 
of  the  sampling.  This  sediment  was  also  sampled  for  size  analysis. 

The  final  measurement  taken  was  the  current  velocity  at  a  point  midway 
between  the  shore  and  the  buoy  line.  The  current  velocity  was  measured 
with  a  Universal  Current  Meter  10.002. 

Due  to  the  location  of  the  sediment  to  be  dredged  at  this  site, 
the  discharge  line  from  the  U.S.  Army  Corps  of  Engineer's  Dredge 
Thompson  to  the  shore  ran  across  the  entire  width  of  the  main  channel. 
This  effectively  blocked  any  barge  traffic  from  either  direction.  For 
this  reason,  the  dredging  operation  was  frequently  halted  and  the  dis¬ 
charge  line  separated  for  barges  to  pass.  This  interrupted  the  flow 
from  the  island  and  sampling  was  halted  until  the  flow  was  resumed. 
Sampling  was  not  restarted  for  a  period  of  time  after  the  flow  had 
returned  to  allow  time  for  the  plume  to  become  re-established. 

The  dredge  spoil  at  the  Rock  Island,  Illinois  site  was  discharged 
directly  to  the  Illinois  shore  of  the  river  in  what  is  known  as  a 
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"beach  nourishment"  type  of  operation,  see  Figure  4-6.  This  was  the 
only  operation  with  side  bank  disposal.  In  this  type  of  operation,  a 
major  percentage  of  the  discharged  sediment  settles  on  the  river  bank 
while  a  small  portion  of  the  sand  and  the  majority  of  the  silt  and 
clay  fractions  are  retained  in  the  water  that  returns  to  the  river. 
These  fractions  make  up  the  plume. 

Shore  markers  and  buoys  were  again  located  in  such  a  manner  as  to 
encompass  the  plume.  The  grid  marked  out  was  430  m.  long  and 
approximately  100  m.  wide.  Distances  were  taken  with  the  geodimeter 
and  sampling  of  the  plume  was  begun.  The  fluorometer  was  still  in¬ 
operative,  it  would  not  hold  a  zero  reading,  so  many  discrete  samples 
were  taken  to  be  analyzed  for  suspended  solids  and  turbidity  at  the 
laboratory.  All  samples  were  taken  at  the  three  foot  depth,  which  was 
approximately  mid-depth.  Water  samples  were  taken  at  several  points 
in  the  plume  for  size  analysis  of  the  suspended  sediments.  Samples 
were  taken  of  the  water  flowing  across  the  bank  before  entering  the 
river,  the  water  at  the  head  of  the  plume,  and  water  approximately 
100  meters  downstream  from  the  head  of  the  plume.  A  sample  of  the 
deposited  sediment  near  the  dredge  discharge  was  also  collected  for 
size  analysis. 

Discharge  flow  at  this  location  was  also  quite  intermittent.  Due 
to  the  morphometry  of  the  river  bed,  there  were  times  when  very  little 
sediment  was  being  dredged  and  discharged.  These  periods  of  pure  water 
discharge  could  last  for  minutes.  During  these  times,  very  little 
suspended  material  was  being  added  to  the  plume.  The  dredging 
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operation  was  also  halted  several  times  to  move  the  discharge  pipe 
further  upstream.  This  had  the  dual  effect  of  stopping  the  sampling 
activity  and  moving  the  source  of  the  plume  to  a  new  location. 

Field  Results 

The  results  from  the  Hannibal,  Missouri  sampling  trip  were  size 
analyses  on  suspended  material  in  two  samples.  The  first  was  a  sample 
of  water  flowing  very  near  to  the  dredge  discharge.  It  was  attempted 

<y 

to  get  a  homogenous  sample  of  material  being  discharged  from  the 
dredge  but  this  was  not  possible.  When  material  was  discharged  from 
the  dredge  discharge  pipe,  a  large  portion  of  the  solids  immediately 
settled.  The  water  portion  of  the  dredged  material  flowed  over  this 
mounded  sand.  This  water  was  sampled  for  size  analysis.  The  results 
of  the  size  analysis  are  shown  in  Figure  B-l,  Appendix  B  and  a 
summary  is  shown  in  Table  4-1.  This  sample  contained  2100  mg/1  of 
suspended  solids. 

The  second  sample  was  of  water  flowing  overland  across  Armstrong 
Island.  This  sample  was  collected  approximately  one  half  mile  from 
the  discharge  point.  The  water  was  fairly  slow  moving  and  had  passed 
through  some  relatively  quiescent  pools.  It  was  felt  that  this  water 
was  indicative  of  the  water  that  would  have  returned  to  the  river,  had 
there  been  return  flow. 

The  suspended  solids  content  of  this  sample  was  74  mg/1.  It  can 
be  seen  in  Table  4-1  that  the  size  of  the  suspended  solids  in  the  over¬ 
land  flow  water  was  much  smaller  than  the  size  of  the  suspended 


Results  of  size  analysis  of  water  and  sediment  samples. 
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solids  in  the  discharged  material.  This  is  due  to  settling  of  the 
larger  material  in  the  quiescent  pools.  The  results  of  the  size  anal¬ 
ysis  can  be  seen  in  Figure  B-2,  Appendix  B. 

The  size  analyses  were  performed  by  Mr.  Wilbur  Matthes,  Jr.,  United 
States  Geological  Survey,  Iowa  City.  The  particle  sizes  were  analyzed 
by  the  Visual  Accumulation  Method  and  the  Pipet  Method.  The  Visual 
Accumulation  method  gives  an  analysis  in  the  range  of  62  to  1000  micro¬ 
meters  (ym).  The  Pipet  Method  gives  an  analysis  of  the  particles  in 
the  range  of  2  to  62  ym.  One  sample  was  analyzed  by  the  Dry  Sieve 
Method  which  gives  an  analysis  of  particles  in  the  range  of  62  to 
4000  ym. 

The  mid-depth  suspended  solids  plume  as  sampled  at  Keithsburg  is 
shown  in  Figure  4-7.  It  can  be  seen  that  the  plume  hugged  the  shoreline 
and  exhibited  little  lateral  dispersion;  the  plume  is  less  than  20 
meters  wide. 

The  results  from  the  four  size  analyses  performed  are  shown  sum¬ 
marized  in  Table  4-1  and  in  Figures  £-3  through  B-6,  Appendix  B.  The 
samples  analyzed  were  island  mud,  deposited  sand,  discharge  creek  and 
head  of  plume.  The  first  two  were  sediment  samples  while  the  second 
two  were  water  samples.  The  island  mud  sample  was  the  material  being 
eroded  to  form  the  suspended  solids  plume.  It  was  mostly  silt  with  a 
small  amount  of  sand  an'd  clay.  The  deposited  sand  is  material  de¬ 
posited  as  the  runoff  water  entered  the  river.  Of  the  two  water 
samples,  the  first  was  taken  in  the  eroded  discharge  creek  before  enter¬ 
ing  the  river,  and  the  second  was  taken  at  the  head  of  the  plume,  after 
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4-7.  Mid-depth  suspended  solids  plume  at  Keithsburg,  Illinois 
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the  sand  had  been  deposited.  This  can  be  seen  by  comparing  the  two 
particle  size  frequency  plots;  the  head  of  plume  sample  shows  a  lower 
percentage  of  sand  than  the  discharge  creek  plot. 

The  mid-depth  suspended  solids  plume  as  sampled  at  Rock  Island 
is  shown  in  Figure  4-8,  with  iso-concentration  lines.  It  can  be  seen 
that  this  plume  also  hugged  the  shoreline  and  exhibited  little  lateral 
dispersion  over  500  m.  downstream  distance.  The  drop  in  suspended 
solids  concentration  between  200  m.  and  350  m.  is  assumed  to  be  caused 
by  a  prolonged  period  of  low  solids  concentration  in  the  discharge. 

The  sampling  was  discontinued  at  450  m.  because  of  a  large  widening 
and  change  in  river  morphometry  at  this  point. 

There  were  four  samples  collected  at  Rock  Island  for  size  analysis, 
one  sediment  sample  and  three  water  samples,  see  Table  4-1  and  Fig¬ 
ures  B-7  through  B-10,  Appendix  B.  The  sediment  sample  was  of  sediment 
near  the  dredge  discharge, but  away  from  the  bank  approximately  5  feet. 
This  sediment  is  material  that  had  been  dredged  from  the  channel, 
discharged  on  the  bank  and  carried  back  into  the  river  by  the 
water.  It  can  be  seen  that  this  was  very  large  material.  A  sample 
was  collected  near  the  discharge  pipe  in  the  same  manner  as  the  sample 
collected  at  the  Hannibal,  Missouri  dredge  operation.  It  can  be 
seen  that  these  samples  are  very  similar.  The  water  samples  collected 
at  the  head  of  the  plume  and  100  meters  downstream  in  the  plume 
show  interesting  results.  The  suspended  solids  in  the  plume  at  the 
head  consist  of  primarily  sand,  while  only  100  meters  downstream, 
there  is  almost  no  sand.  Another  interesting  observation  is  that  the 


41 


LATERAL  DISTANCE  (  meters  ) 


Figure  4-8.  Mid-depth  suspended  solids  plume  at  Rock  Island 
Illinois  site. 
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silt  and  clay  fractions  have  reversed,  there  being  a  much  higher  per¬ 
centage  clay  in  the  body  of  the  plume  than  at  the  head  of  the  plume. 

The  turbidity  of  the  samples  collected  for  size  analysis  was  mea¬ 
sured  ir.  the  laboratory  with  a  Hach  Model  2100  Turbidimeter.  The  max¬ 
imum  turbidity  measured  was  33  nephelometric  turbidity  units  (NTU) 
above  ambient  at  the  plume  source.  The  ambient  turbidity  was  22  NTU. 
The  turbidity  in  the  plume  rapidly  decreased  with  downstream  distance; 
the  turbidity  had  decreased  to  15  NTU  at  100  m.  downstream.  Figure  0-2 
in  Appendix  D  shows  the  relationship  between  suspended  solids  and  tur¬ 
bidity  for  the  Rock  Island  samples.  The  correlation  coefficient  for 
this  data  is  0.87. 
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CHAPTER  V 

MODEL  RESULTS  AND  DISCUSSION 


Graphical  Solution 

The  model  development  by  Schubel ,  Carter  et  al.,  (197S)  was  first 
used  to  try  to  simulate  the'  observed  field  data.  It  was  decided  to 
simulate  the  data  from  the  Rock  Island,  Illinois  sampling  trip. 

Model  Input  Parameters 

There  are  six  input  parameters  to  the  model;  a)  the  rate  of 
addition  of  suspended  solids  to  the  receiving  water,  b)  the  average 
vertical  thickness  of  the  plume,  c)  the  mean  particle  settling  veloc¬ 
ity,  d)  the  diffusion  velocity,  e)  the  time  interval  for  the  plume 
to  reach  its  maximum  length  and  f)  the  average  current  velocity  of 
the  receiving  water.  Each  of  these  parameters  will  be  discussed  as 
pertaining  to  the  Rock  Island,  Illinois  site. 

Rate  of  addition  of  suspended  solids  to  the  receiving  water  (q) 

The  rate  of  addition  of  suspended  solids  to  the  receiving  water 
is  a  function  of  the  size  of  the  dredge,  the  type  of  material  being 
dredged  and  the  amount  of  time  for  settling  before  the  discharged 
water  returns  to  the  receiving  water.  Since  the  operation  at  Rock 
Island  was  side  bank  disposal,  there  was  essentially  no  time  for 
settling  before  the  discharged  water  re-entered  the  river.  The  amount 


of  suspended  material  entering  the  river  and  the  rate  of  addition  can 
be  calculated  in  several  ways. 

The  fraction  of  the  total  solids  discharged  from  the  dredge  that 
becomes  incorporated  into  the  plume  has  been  calculated  to  vary  from 
1%  to  5%  (Schubel ,  Carter  et  a!.,  1978).  The  mass  of  material  dis¬ 
charged  from  the  dredge  per  unit  time,  Q  ,  can  be  calculated.  At 
Rock  Island,  11,596  cubic  yards  of  material  were  dredged  in  an  operat¬ 
ing  time  of  18.58  hours  (personal  communication  with  Mr.  Dick  Baker, 
Chief  of  Operations,  Rock  Island  District,  U.S.  Army  Corps  of  Engineers). 
Using  these  values  and  assuming  the  sediment  to  be  85%  solids, 

O 

Qm  =  2.39  x  10  mg/sec.  The  fraction  remaining  suspended  and  becoming 
incorporated  into  the  plume  is  assumed  to  be  the  silt  and  clay  fraction, 
which  from  Table  4-1,  is  seen  to  be  5.0%  at  the  beginning  of  the  plume. 
Therefore, -the  rate  of  addition  of  suspended  particulates  to  the  plume, 
q,  is  equal  to  1.20  x  107  mg/sec. 

An  alternate  method  of  calculation  of  the  rate  of  addition  of 
suspended  solids  to  the  plume  is  to  calculate  the  value  of  q  =  ux  AC^, 
where  A  is  the  cross-sectional  area  of  the  head  of  the  plume,  C^  is  the 
concentration  of  suspended  solids  at  the  head  of  the  plume  and  ux  is 
the  mean  plume  velocity  in  the  longitudinal  direction.  From  Figure  4-8, 
it  is  seen  that  the  width  of  the  plume  is  approximately  50  meters  at 
the  source,  and  the  concentration  is  approximately  112  mg/1  at  that 
point.  The  average  depth  of  the  river  was  measured  to  be  6  feet  and 
the  mean  current  velocity  was  0.40  meters/second.  Using  this 


information: 

A  =  (50  m) (2  m)  =  100  m2  (5,1) 

q  =  (0.40  m/sec) (100  m2)(112  mg/l)(1000  1/m3)  =  4.48  x  106  mg/sec. 

(5.2) 

It  is  seen  that  there  is  a  large  disagreement  in  q  calculated  by 
the  two  methods.  Since  the  objective  is  to  try  to  match  the  observed 
suspended  solids  plume,  the  value  q  =  4.5  x  10^  mg/sec  is  chosen. 
Evidently  some  of  the  silt  and  clay  must  settle-out  in  a  dense  wedge 
as  the  discharge  water  first  enters  the  river.  Approximately  2%  of 
the  total  sediment  that  is  dredged  actually  enters  the  River  and  be¬ 
comes  entrained  in  the  plume. 

Average  vertical  thickness  of  the  plume,  D 

The  depth  of  the  river  was  measured  at  several  locations  in  the 
suspended  solids  plume.  The  average  depth  was  determined  to  be  approx¬ 
imately  6  feet.  Schubel,  Carter  et  al . ,  (1978)  advise  using  a  value 
of  one  half  the  total  water  depth  in  areas  where  the  water  depth  is 
8  feet  or  less.  Therefore,  the  value  D  =  3  feet  =  0.9  meters  is  chosen. 

Mean  particle  settling  velocity,  W 

The  mean  particle  size  can  be  determined  from  the  size  analysis 
on  the  suspended  solids.  Since  the  sand  settles  immediately,  the 
material  forming  the  plume  is  the  silt  and  clay  fraction.  The  mean 
particle  size  of  the  silt/clay  fraction  was  determined  to  be  0.02  mm. 
Using  Stoke's  Law  and  a  water  temperature  of  50°  F,  the  mean  particle 


settling  velocity  was  calculated,  W  =  0.027  cm/sec. 
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Diffusion  velocity,  oo 

Schubel,  Carter  et  al.,  (1978)  reported  the  range  of  the  longi¬ 
tudinal  and  lateral  diffusion  velocity  in  open  rivers  to  be  0.2  - 
0.5  cm/sec.  The  value  of  0.5  cm/sec  was  chosen. 

Time  interval  for  the  plume  to  reach  its  maximum  length,  t 

The  maximum  length  of  a  suspended  solids  plume  in  a  river  is 
determined  by  the  settling  velocity  of  the  suspended  particle  and  the 
vertical  distance  the  mean  particle  must  settle  (Barnard,  1978).  For 
the  Rock  Island  case,  W  =  0.027  cm/sec  and  D  =  3  feet, 

‘  •  ff  *  '  3370  seconds  (5-3) 


Average  current  velocity  of  the  receiving  water,  u.. 

The  current  velocity  of  the  river  was  measured  at  several  locations 
within  the  suspended  solids  plume.  The  average  current  velocity  was 
calculated  to  be,  ux  =  0.4  m/sec. 

Non-dimensional  Ratios  and  Scaling  Factors 
The  suspended  solids  model  presented  by  Schubel,  Carter  et  al., 
(1978)  is  in  the  form  of  a  series  of  graphs;  The  graphs  were  developed 
as  functions  of  the  following  non-dimensional  ratios  and  scaling 
factors. 


Ratio  of  diffusion  velocity  to  advective  velocity,  ^/u 

_  0.5  cm/sec  _  n 
ux  '  40  cm/sec  "  U'UI'3 


(5.4) 
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This  ratio  indicates  that  the  longitudinal  dispersion  is  small  in 
comparison  to  the  mean  longitudinal  velocity. 

Ratio  of  the  plume  age  to  the  settling  time,  y 

Wt  (0.027  cm/sec) (3370  sec)  _  ,  tc  c\ 

Y  “  D  "  91  cm  1 

The  value  of  ywill  always  be  equal  to  1  in  a  river  since  t  is  de¬ 
fined  as  D/W. 


Distance  Scaling  Factor,  DSF 

DSF  =  uxt  =  (0.4  m/sec)(3370  sec)  =  1350  m  (5-6) 

This  is  the  expected  distance  of  travel  for  the  mean  particle 
which  falls  from  the  surface  to  the  bottom. 


Concentration  Scaling  Factor,  CSF 


CSF 


q  _  2.24  x  106  mg/sec  (1000  cm3/!) 
Trcobt  (it)  (o!  5  cm/sec) *  (91  cm)  (3370  sec) 


=  18,600  mg/1  (5.7) 


Calculation  of  the  Centerline  Concentrations 
The  above  ratios  and  factors  are  used  along  with  the  graphs  of 
Schubel,  Carter  et  al.,  (1978)  to  calculate  the  concentration  of  the 
suspended  solids  plume  along  the  centerline;  The  centerline  for  a 
sidebank  disposal  operation  in  a  river  is  along  the  near  bank.  The 
model  was  originally  developed  for  estuarine  open  water  disposal,  and 
therefore, no  effects  of  sidebanks  were  included  in  the  solution.  This 
is  easily  modified  for  sidebank  disposal  in  a  river  by  assuming  the 
bank  is  a  reflecting  barrier.  The  effect  of  this  reflecting  barrier 
on  a  plume  resulting  from  sidebank  disposal  can  be  described  as  folding 
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the  plume  back  on  itself  along  the  ceterline.  The  net  effect  is  that 
the  suspended  solids  concentrations  calculated  with  this  model  must  be 
doubled  to  describe  sidebank  disposal  in  a  river. 

Suspended  solids  concentration  at  distance  u^t 

The  first  step  in  determining  the  suspended  solids  concentrations 
along  the  centerline  is  to  determine  the  suspended  solids  concentration 
at  distance  uxt.  This  concentration  is  found  by  using  Figure  5-1, 
(Barnard, 1978) .  Enter  Figure  5-1  at  the  calculated  value  of  w/u  . 

A 

Move  vertically  to  the  curve  corresponding  to  the  calculated  value  of 

Y  and  horizontally  to  determine  the  value  of 

Concentration,  mq/1  at  distance  u..t  =  0.0045 
CSF  * 

Therefore, the  suspended  solids  concentration  at  1350  m  is  calculated 
to  be  equal-  to  84  mg/1  above  the  ambient  river  value.  Doubling  this 
value  to  account  for  reflection  from  the  bank  gives  a  value  of  167  mg/1 
above  ambient. 

Distance,  x,  where  centerline  concentration  is  a  specified  concentration 
above  ambient 

The  next  step  in  determining  the  suspended  solids  concentrations 
along  the  centerline  is  to  choose  a  centerline  concentration  and  find 
the  distance  downstream  that  corresponds  to  this  concentration.  As  an 
example,  the  distance  where  the  centerline  suspended  solids  concentra¬ 
tion  =  1000  mg/1  above  ambient  will  be  calculated. 


CONCENTRATION,  mq/4  AT  DISTANCE  U3 
"  CSF  ~ 
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Figure  5-1.  Relationship  between  and  Solids  concentration  at  _di s t-snce  X 
for  y  =  0.1,  1 ,  3.2  and  10. 
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,  .  .  _ 1000  mg/1 _ _  1000  mq/1  _  ,  n  ,  c  0\ 

1)  Calcul  te.  concentration  at  distance  uxt  167  mg/1  '  ' 

2)  Use  this  ratio  to  enter  Figure  5-2  from  Barnard  (1978)  along 

the  ordinate.  Move  horizontally  to  the  curve  corresponding 

to  a)/ux  and  then  vertically  to  determine  the  value  of 

Distance  x 
DSF 


Figure  5-2  is  for  Y  =  1.  Figures  for  y  =  0.01,  0.1,  10  and 
100  are  included  in  Schubel,  Carter  et  al.,  (1978).  Multiply¬ 
ing  this  value  by  DSF  gives  the  distance  at  which  the  center- 
line  suspended  solids  concentration  is  1000  mg/1  above 


ambient. 


For  the  Rock  Island  site,  co/u  =  0.013  and  DSF  =  1350  m. 

A 

Figure  5-2  shows 

-  Distance^  =  0-52  (5.9) 


Therefore,  the  distance  where  the  centerline  suspended  solids 
concentration  is  1000  mg/1  above  ambient  is  equal  to  700  m. 

3)  Steps  1  and  2  are  repeated  for  as  many  different  suspended 
solids  concentrations  as  are  needed  to  adequately  describe 
the  centerline  of  the  plume.  Values  calculated  for  the 
Rock  Island  site  are  shown  in  Table  5-1. 

It  can  be  seen  in  Figure  5-2  that  the  curve  for  co/ux  is  nearly 
vertical  below 


Concentration  at  distance  x  _  -j 

Concentration  at  distance  uxt 

For  this  reason,  plume  concentrations  can  not  be  calculated  at  distances 
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Concentration 
(mg/1 ) 

Concentration 

Distance  x 

Distance  x 
(m) 

Concentration  at  distance  u  t 

A 

i 

DSF 

25,000 

150 

0.032 

43 

10,000 

60 

0.085 

115 

5,000 

30 

0.14 

189 

2,500 

15 

0.23 

310 

1 ,000 

6 

0.52 

500 

3 

0.80 

1080 

167 

1  1 

1.00 

1350 

Table  5-1.  Downstream  distance  corresponding  to  various  suspended 
solids  concentrations  along  the  centerline.  Rock  Island,  Illinois  site. 
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beyond  uxt.  For  the  Rock  Island  plume,  the  plume  can  not  be  described 
at  a  distance  beyond  1350  m.  The  suspended  solids  concentration  at 
this  point  is  167  mg/1. 


Lateral  Dimensions  of  the  Plume 

The  plume  described  by  this  model  is  approximately  Gaussian  and 
therefore,  the  lateral  dimensions  are  directly  related  to  x,  the  down¬ 
stream  distance.  The  width  of  the  plume  as  determined  by  the  C ( I ) 
isopleth  and  measured  from  the  centerline,  y,  is  determined  by: 


where  a2(x/DSF)  is  determined  from  Figure  5-3,  a2(l)  is  determined 
from  Figure  5-4,  C(x)  =  suspended  solids  concentration  on  the  center- 
line  at  distance  x,  C(I)  =  suspended  solids  concentration  of  the 
isopleth  chosen  to  define  the  plume,  w/ux,  DSF,  CSF  and  x  as  defined 
previously. 

Assume  the  plume  is  defined  by  the  50  mg/1  above  ambient  isopleth. 
For  the  Rock  Island  site,  C(I)  =  25  mg/1  since  the  plume  is  reflected 
from  the  shoreline.  The  width  of  the  plume  can  be  calculated  at  each 
of  the  distances  where  the  centerline  suspended  solids  concentration 
is  known.  To  finish  the  example  calculation,  the  width  of  the  plume 
is  calculated  at  x  =  700  m  or  where  C(x)  =  1000  mg/1  above  ambient. 

(1)  Calculate: 


x  _  700  m 
DSF  1350  m" 


(5.9) 


This  is  the  value  that  was  found  in  Step  2, previously. 


Figure  5-4.  Relationship  between  o2(l)  and  w/u  for  determi 
lateral  dimensions  of  the  plume. 
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(2)  Calculate: 

Ciil  -  25  m-3/J -  =  g  gg  I  o 

CSF  18600  mg/1  u‘uuu 


(5.  11) 


(3)  Use  the  value  x/DSF  to  enter  Figure  5-3  along  the  abscissa. 
Move  vertically  to  the  correct  Vux  curve  and  then  horizon¬ 
tally  to  determine  the  value  of  cr  (x/DSF).  From  Figure  5-3, 
with  x/DSF  =  0.52  and  w/u  =  0.013,  02  (x/DSF)  =  0.23. 

(4)  Use  the  value  w/ux  to  enter  Figure  5-4  along  the  abscissa. 

Move  vertically  to  the  correct  y  curve  and  then  horizontally 
to  determine  the  value  of  cr2(l).  From  Figure  5-4,  with 
oj/ux  =  0.013  and  y  =  1,  o2(l)  =  1. 

(5)  Calculate: 


-X-  s  a 
DSF  \ 


y_x  \  a2 ( i )  /w_\2  r- 

vdsf;  U)  l 


In  C(  I  )/CSF" 
C  x)  J 


(5.10) 


_  (0. 25) (1 ) (0.013)2  (-  in  0.0013\ 

nI  \  1000  ) 

(6)  Calculate: 

y  =  (y/DSF) (DSF)  = (0.024) ( 1350  m)  =  32  m 


=  0.024 


(5.12) 


(7)  Steps  1-6  are  repeated  for  other  values  of  x  and  C(x)  until 
the  shape  of  the  50  mg/1  isopleth’is  adequately  determined. 
Table  5-2  shows  the  values  calculated  for  the  50  mg/1  isopleth 
for  the  Rock  Island  site.  If  other  isopleths  are  desired, 
the  procedure  is  repeated  for  a  different  C(I)  value. 
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Discussion  of  Model 

It  can  be  seen  by  comparing  Table  5-1  to  Figure  4-8  that  the  model 
does  not  predict  the  plume  observed  during  the  dredge  disposal  opera¬ 
tion.  The  maximum  suspended  solids  concentration  observed  was  124  mg/1 
above  ambient  with  the  maximum  extent  of  the  50  mg/1  isopleth  being 
approximately  500  m.  downstream.  The  model  shows  a  maximum  suspended 
solids  concentration  in  excess  of  25,000  mg/1  and  the  maximum  longitud¬ 
inal  extent  of  the  50  mg/1  isopleth  is  greater  than  1350  m. 

Much  of  this  problem  can  be  traced  to  the  assumptions  concerning 
the  type  of  source  of  the  plume.  The  model  assumes  a  point  source 
discharge  which  is  consistent  with  the  mode  of  discharge  in  estuarine 
pipeline  disposal  operations,  but  it  is  not  representative  of  the  side- 
bank  disposal  operations  performed  on  the  Mississippi  River,  GREAT  II 
reach. 

It  can  be  seen  from  Figure  4-8  that  the  concentration  of  suspended 
solids  at  x  =  0  is  approximately  constant  for  a  distance  of  50  m. 

This  line  source  means  that  the  same  amount  of  solids  is  suspended  in 
a  much  greater  volume  of  water  for  the  observed  plume,  as  opposed  to 
the  model  calculated  plume. 

The  solution  to  this  problem  would  be  to  modify  the  model  so  the 
initial  source  condition  would  be  a  line  instead  of  a  point.  Unfortun¬ 
ately,  this  is  not  an  easy  task. 

An  alternative  solution  would  be  to  solve  the  model  assuming  sev¬ 
eral  point  sources  located  at  several  points  across  the  observed  plume 
source  width.  This  type  of  solution  may  have  more  closely  described 
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the  observed  plume,  but  the  problem  of  reflection  from  the  shoreline 
would  have  made  the  solution  very  unwieldy.  It  was  decided  that 
this  type  of  solution  was  beyond  the  scope  of  "a  simple  model"  and  so 
was  unsatisfactory  for  this  study. 

Because  the  model  could  not  be  used  to  predict  the  observed  plume 
at  Rock  Island,  this  model  was  not  used  to  try  to  predict  the  Keithsburg 
plume. 


Analytical  Solution 

The  model  discussed  in  Chapter  III  was  developed  as  an  alternative 
to  the  model  developed  by  Schubel,  Carter,  et  al.  This  analytical 
model  was  developed  to  describe  transport  and  dispersion  of  suspended 
solids  in  a  river.  It  was  decided  to  simulate  the  suspended  solids 
plume  that-was  observed  during  the  Rock  Island  dredging  operation. 

Model  Input  Parameters 

There  are  six  input  parameters  to  the  model,  a)  the  width  of  the 
plume  source,  b)  the  mean  depth  of  the  portion  of  the  river  containing 
the  plume,  c)  the  mean  velocity  of  the  river  in  the  area  of  the  plume, 
d)  the  dispersion  coefficient,  e)  the  settling  velocity  of  the  suspended 
particle  and  f)  the  downstream  distance  from  the  plume  source.  Each 
parameter  will  be  discussed. 

Width  of  the  plume  source,  b 

The  width  of  the  plume  source  is  a  function  of  the  velocity  of 
the  returning  flow,  the  direction  of  that  flow  with  respect  to  the 
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direction  of  the  receiving  water  flow  and  the  velocity  of  the  receiv¬ 
ing  water.  A  plume  resulting  from  sidebank  disposal  would  be  expected 
to  have  a  wider  plume  source  than  a  plume  resulting  from  disposal  at 
a  site  where  the  water  could  not  immediately  return  to  the  river.  A 
plume  developing  in  a  backwater  area  characterized  by  a  slow  moving 
current  would  be  expected  to  have  a  wider  plume  source  than  a  plume 
developing  near  the  main  channel  of  a  river  with  a  fast  current. 

Mean  depth  of  the  plume,  D 

The  plumes  resulting  from  dredge  disposal  operations  are  generally 
shore-attached  and  are  therefore  in  areas  of  varying  depth.  The  model 
assumes  constant  depth.  Therefore,  an  attempt  should  be  made  to  mea¬ 
sure  the  depth  at  several  locations  in  the  area  of  the  plume  to  deter¬ 
mine  an  average  depth.  This  average  depth  will  be  used  in  the  model, 
di rectly. 

Mean  river  velocity,  u^ 

The  mean  velocity  of  the  river  in  the  immediate  vicinity  of  the 
plume  must  be  known. 

Dispersion  coefficient, 

The  lateral  dispersion  coefficient  must  be  either  measured  or  cal¬ 
culated  from  empirical  relationships.  Assuming  Elder's  (1959)  relation¬ 
ship  of  Ky  =  0.23  DIJ*,  Ky  was  calculated  to  be  approximately  100  cm2/sec. 
Allowing  for  some  effects  of  the  sloping  channel  bottom,  the  value  of 
^  was  chosen  to  be  300  cm2/sec. 
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Settling  velocity  of  the  suspended  particle,  W 

The  terminal  settling  velocity  of  a  suspended  particle  is  deter¬ 
mined  by  its  specific  gravity  and  its  size.  The  particle  size  can  be 
estimated  by  fall  velocity  analysis.  Knowing  the  particle  size  and 
assuming  a  specific  gravity,  Stoke's  Law  can  be  used  to  determine  the 
terminal  settling  velocity.  A  chart  is  included  in  Barnard  (1978)  and 
is  reproduced  in  Appendix  D,  relating  particle  size  to  terminal  settling 
velocity.  This  chart  can  be  used  for  settling  velocity  approximations. 

Downstream  distance  from  the  plume  source,  x 

The  model  equation  given  in  Chapter  III  is  solved  at  a  particular 
distance  downstream  from  the  plume  source.  To  determine  the  entire 
plume,  the  model  must  be  solved  several  times  with  different  x  distances. 

Solution  to  Rock  Island,  Illinois  Plume 
A  solution  was  first  attempted  using  a  mean  settling  velocity,  as 
with  the  model  of  Schubel ,  Carter,  et  al.,  (1978).  It  was  not  possible 
to  calculate  a  plume  that  resembled  the  field-observed  plume  using  this 
technique.  Therefore  it  was  decided  to  calculate  the  plume  resulting 
from  each  of  three  different  size  fractions  and  sum  the  individual 
concentrations  to  get  the  overall  plume. 

The  final  solution  involved  three  size  fractions;  sand,  silt  and 
clay.  The  proportion  of  each  fraction  was  determined  by  knowing  the 
size  analysis  at  a  point  in  the  plume  and  the  concentration  of  total 
solids  in  the  plume  at  various  locations.  Table  4-1  shows  that  the 
silt/clay  ratio  was  approximately  25:75  at  a  point  in  the  plume. 
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Figure  4-8  shows  a  suspended  solids  concentration  of  approximately 
45  mg/1  at  the  end  of  the  plume.  Assuming  no  clay  settled  out  of  the 
plume  over  the  short  length  of  the  measured  plume,  the  initial  clay 
concentration  was  calculated  to  be  approximately  35  mg/1  or  30%  of  the 
initial  suspended  solids  concentration.  The  proportion  of  sand  was  de¬ 
termined  from  Figure  4-8.  The  initial  suspended  solids  concentration 
was  approximately  112  mg/1  and  the  suspended  solids  concentration  at 
100  m  was  approximately  60  mg/1.  Assuming  all  of  the  sand  had  settled 
in  the  first  100  m  and  that  little  of  the  silt  and  clay  had  been  re¬ 
moved,  gave  an  initial  sand  concentration  of  52  mg/1  or  approximately 
45%  of  the  initial  suspended  solids.  The  remaining  25%  of  the  initial 
suspended  solids  was  assumed  to  be  the  silt  fraction.  Thus,  the  com¬ 
position  of  the  suspended  solids  at  the  plume  source  was  approximately 
45%  sand,  25%  silt  and  30%  clay. 

The  settling  velocity  for  each  of  these  fractions  was  estimated 
from  Appendix  D  and  from  the  size  analysis  of  the  material  entering 
the  plume, Figure  B-9,  Appendix  B.  The  mean  diameter  of  the  sand 
fraction  was  determined  to  be  0.26  mm,  corresponding  to  a  settling  vel¬ 
ocity  of  0.02  m/sec.  The  silt  fraction  mean  particle  size  was  0.026  mm 
with  a  settling  velocity  of  0.003  m/sec.  The  settling  velocity  of  the 
clay  fraction  was  chosen  as  0.000001  m/s. 

The  following  parameters  were  used  as  input  to  solve  the  model 
for  the  Rock  Island  simulation: 

b  =  25  m 

D  =  2  m 
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ux  =  0.4  m/s 
Ky  =  0.03  m2/s 

Wsand  =  °-02  m/s 
Wsilt  =  m/s 

Wclay  =  °* 000001  m/s 
Q  C 

-R-2-  =  112  mg/1 


where  QqCo  is  the  concentration  of  suspended  material  at  x  =  0. 


These  values  are  used  in  Equation  3.24  to  simulate  the  plume. 


Qb  C  (y,x)  _ 

F 

Lb\  F  Ar 

exp 

Wx 

Q  C 
v0  0 

V  c 

’y)  V 

°y 

Du 

x_ 

(3.24) 


The  method  of  solution  to  Equation  3.24  follows. 

(1)  -Choose  a  distance,  x,  downstream  from  the  source.  For 
example,  choose  50  m. 

(2)  Calculate  the  value  of  the  exponential  term. 


exp 


[- 


Wx  1  -  pvn  r  (0.02  rn/s) (50  m) 
Dux  '  exp  '  - 

L  * . 


(2  m) (0.4  m/s) 


=  0.287 


(3)  Choose  value  of  y/o  .  This  value  corresponds  to  the  distance 
from  the  plume  centerline,  y.  The  plume  centerline,  y  =  0, 
is  defined  as  the  shoreline  along  which  the  plume  develops. 

To  calculate'  the  centerline  concentration,  y/o^  =  0. 

(4)  Calculate  the  value  of  o  . 
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Appendix  C. 

F  (9.09)  =  1.0 
F  (-9.09)  =  0.0 

(10)  Insert  values  calculated  in  steps  2  and  10  into  Equation  3.24 
to  calculate  the  proportion  of  chosen  size  fraction  remaining 
in  the  plume  at  point  (y,x). 


Qb  C  (y,x) -=  (1.0  -0.0)  (0.287)  =0.287. 

Q  C 
yo  o 
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(11)  Calculate  the  initial  concentration  of  the  sediment  fraction 
chosen. 

%  sediment  chosenj  =  (112  mg/1) (0.45  sand) 

=50.4  mg/1  sand 

(12)  Calculate  the  concentration  of  the  chosen  sediment  fraction 
at  point  (y ,x) . 

C  (y,x)  =  (50.4  mg/1 ) (0. 287)  =  14.5  mg/1 

(13)  Repeat  steps  3  through  12  for  a  sufficient  number  of  values 
of  y /a  to  determine  the  concentration  of  the  chosen  sedi¬ 
ment  fraction  in  the  plume  cross-section  at  distance  x. 

(14)  Repeat  steps  2  through  13  for  the  various  sediment  fractions. 

(15)  Sum  the  values  calculated  in  step  12  for  each  point  (y,x)  to 
determine  the  overall  suspended  solids  concentration  at 
point  (y,x). 

(16)  Repeat  steps  1  through  15  for  sufficient  number  of  values  of 
x  to  determine  the  dimensions  of  the  plume  and  the  concen¬ 
trations  in  the  plume. 

The  simulated  suspended  solids  plume  calculated  for  the  Rock 
Island,  Illinois  site  is  shown  in  Table  5-3. 

Figures  5-5,  5-6  and  5-7  all  show  the  simulated  suspended  solids 
plume.  Figure  5-5  shows  the  plume  superimposed  on  the  field  data. 
Figures  5-6  and  5-7  show  the  simulated  plume  to  a  distance  of  10,000 
meters  for  different  values  of  Ky.  Figure  5-6  is  for  Ky  =  300  cm2/sec., 
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5-7.  Far  field  suspended  solids  plume  for  Rock  Island,  Illinois 
site.  Ky  =  1000  cmvsec. 
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the  same  value  used  for  Figure  5-5,  while  Figure  5-7  is  for 
Ky  =  1000  cm2/sec.  Since  the  far  field  plume  was  not  measured  during 
the  field  sampling,  there  is  no  basis  for  choosing  either  one  as  the 
correct  plume.  However  the  turbulence  scale  would  increase  as  the 
lateral  dimensions  of  the  plume  increases  downstream,  so  one  might 
expect  an  increase  in  K  ,  the  lateral  dispersion  coefficient. 

Solution  to  Keithsburg,  Illinois  Plume 
The  observed  suspended  solids  plume  at  Keithsburg,  Illinois  was 
quite  different  from  the  plume  observed  at  Rock  Island,  Illinois.  It 
can  be  seen  from  Figure  4-7  that  the  source  width  is  much  smaller;  it 
was  estimated  to  be  3  meters.  Table  4-1  indicates  that  a  very  dif¬ 
ferent  sand:silt:clay  ratio  was  measured  at  Keithsburg;  approximately 
2:67:31  at  the  head  of  the  plume.  Figure  B-6  in  Appendix  B  shows 
mean  diamters  for  the  particles  to  be;  sand  =  0.0086  cm  and 
silt  =  0.0017  cm.  These  correspond  to  settling  velocities  of  0.005  m/sec 
for  sand  and  0.00022  m/sec  for  silt.  The  settling  velocity  for  clay 
was  again  chosen  to  be  1  x  10"6  m/sec.  The  stream  velocity  was  mea¬ 
sured  to  be  0.35  m/sec.  The  initial  suspended  solids  concentration 
was  estimated  to  be  75  mg/1  from  Figure  4-7.  The  parameters  used  as 
input  to  the  model  for  the  Keithsburg,  Illinois  simulation  are  summar¬ 
ized  below. 

b  =  3  m 
0  =  2m 
ux  =  0.35  m/s 
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Ky  =  0.03  m2/s 

Wsand  =  °*005  m/s 

Wsilt  =  °* 00022  m/s 

Wc1ay  =  °* 000001  m/s 
Q  C 

-s-2-  =  75  mg/1 


The  simulated  suspended  solids  plume  for  Keithsburg,  Illinois  is  shown 
in  Table  5-4.  Figure  5-8  shows  this  model  plume  superimposed  on  the 
field  data  while  Figure  5-9  shows  this  same  plume  in  the  far  field. 
Figure  5-10  shows  what  the  far-field  plume  might  look  like  if 
Ky  =  1000  cnr/sec.  Again,  the  far-field  plume  was  not  measured  in  the 
field  so  it  can  not  be  determined  which  plume.  Figure  5-9  or  Figure 
5-10  is  more  correct. 


Discussion  of  the  Analytical  Model 
Inspection  of  the  plumes  generated  by  the  model  (Figures  5-5 
through  5-10)  yield  several  physical  parameters  that  can  be  estimated 
by  solving  the  model  for  a  particular  dredge  disposal  operation.  Some 
of  these  parameters  are,  the  plume  centerline  suspended  solids  concen¬ 
tration,  the  lateral  suspended  solids  concentrations,  the  amount  of 
solids  being  deposited  at  some  point  in  the  receiving  river,  the  di¬ 
lution  volume  for  dissolved  substances,  and  the  maximum  length  of 
the  plume. 

The  model  does  a  good  job  of  predicting  the  centerline  (near  bank) 
suspended  solids  concentration.  Figure  5-11  and  Figure  5-12  show  the 


Table  5-4  (continued) 
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Figure  5-8.  Simulated  suspended  solids  plume  for  Keithsburg,  Illinois 
site  superimposed  on  field  data. 
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Figure  5-9.  Far  field  suspended  solids  plume  for  Keithsburg,  Illinois 
site.  K.  =  300  cm2/sec. 
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Figure  5-11.  Model  prediction  vs.  field  observation  for  centerline 
Rock  Island  site. 
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suspended  solids  concentration  predicted  and  observed  as  a  function  of 
downstream  distance  for  the  Rock  Island  and  Keithsburg  sites,  respect¬ 
ively.  Figure  5-11  shows  good  agreement  between  the  predicted  and 
observed  suspended  solids  concentrations  for  the  first  150  meters  down¬ 
stream.  The  low  field  measurements  observed  between  150  and  350  meters 
are  assumed  to  be  due  to  an  extended  period  of  dredging  in  an  area 
of  deep  water,  hence,  lower  than  normal  solids  concentration  being 
discharged  from  the  dredge.  At  the  end  of  the  observed  plume,  there 
is  about  a  5  mg/1  difference  between  the  observed  concentration  and 
the  predicted  concentration.  Figure  5-12  shows  excellent  agreement 
between  the  observed  and  the  predicted  suspended  solids  concentration 
at  Keithsburg. 

The  model  is  less  successful  in  predicting  the  degree  of  lateral 
dispersion.*  It  can  be  seen  in  Figure  5-5  for  Rock  Island  that  consid¬ 
erable  suspended  solids  concentrations  were  observed  beyond  the  10  mg/1 
and  1  mg/1  isopleths  predicted  by  the  model.  It  can  also  be  seen  that 
there  were  few  samples  taken  inside  the  area  bounded  by  the  10  mg/1 
predicted  isopleth,and  that  those  samples  that  were  taken  show  a  great 
deal  of  variability  with  no  smooth  concentration  gradient. 

In  Figure  5-8  for  Keithsburg,  no  samples  were  taken  within  the 
predicted  plume  at  any  point  other  than  the  centerline.  Those  samples 
taken  beyond  the  predicted  plume  show  no  excess  suspended  solids  in  the 
stream  resulting  from  dredge  disposal. 

In  order  to  assess  the  effect  of  dredge  spoils  disposal  on  the 
benthic  community,  the  amount  of  dredged  material  deposited  on  the 


bottom  must  be  calculated.  This  can  be  readily  calculated  with  the 
known  parameters  and  the  predicted  suspended  solids  plume. 

At  steady  state,  the  rate  of  solids  deposition  at  any  point  in  the 
plume  is: 

Deposition  Rate  =  C  (y,x)  KjD  (5.13) 

where  C  (y,x)  is  the  mean  suspended  solids  concentration  at  the  point 
(y,x),  K$  is  the  settling  rate  constant  and  D  is  the  depth.  The 
settling  rate  constant  can  be  calculated  by  w/q,  where  W  is  the  settling 
velocity  and  D  is  the  depth.  The  solids  deposition  rate  can  now  be 
expressed: 

Deposition  Rate  =  C  (y,x)  W  (5.14) 

Equation  5.14  is  solved  for  each  sediment  fraction  and  the  results  are 
summed  for  the  total  solids  deposition  rate.  The  concentration  of 
each  sediment  fraction  at  many  points  in  the  plume  is  given  in  Table  5-3 
for  Rock  Island  and  Table  5-4  for  Keithsburg.  Multiplying  these  concen¬ 
trations  by  their  respective  settling  velocities  gives  the  deposition 
rate  for  that  sediment  fraction.  The  results  of  this  calculation  are 
given  in  Table  5-5  and  Figure  5-13  for  Rock  Island  and  Table  5-6  and 
Figure  5-14  for  Keithsburg. 

To  determine  the  total  mass  of  solids  deposited  at  any  point, 
multiply  the  deposition  rate  at  that  point  by  the  time  of  operation 
of  the  dredge.  The  depth  of  solids  deposited  at  any  point  can  then  be 
calculated  by  assuming  a  solids  density  and  a  percent  solids.  The 
depth  of  solids  deposited  at  all  points  in  the  plume  at  Rock  Island  is 
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Table  5-6  (continued). 


LATERAL  DISTANCE,  Y  (meters) 


Figure  5-14.  Deposition  rate  at  all  points  in  plume  for 
Keithsburg,  Illinois. 
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shown  in  Figure  5-15.  The  solids  density  is  assumed  to  be  2.65  gm/cm3 
and  the  sediment  is  assumed  to  be  85%  solids.  It  can  be  seen  in  Fig¬ 
ure  5-15  that  only  at  the  head  of  the  plume  is  there  significant  solids 

deposition.  At  a  distance  of  approximately  50  m  downstream,  only 
1  mm  of  sediment  accumulates  during  the  18.58  hours  of  dredging.  It 
should  be  noted  that  this  calculation  does  not  include  sediment  which 
immediately  falls  to  the  bottom  in  a  dense  wedge  near  the  head  of  tlfe 
plume. 

The  depth  of  solids  deposited  from  the  plume  at  Keithsburg  is 
shown  in  Figure  5-16.  It  can  be  seen  that  there  is  less  than  0.5  mm  of 

sediment  deposited  in  the  16  hours  of  dredging. 

Another  area  of  interest  in  impact  assessment  of  dredge  spoils 
discharge  operations  is  desorption  of  substances  previously  adsorbed 
to  sediment  particles  during  dredging  and  disposal.  After  desorption, 
these  substances  are  dissolved  in  the  discharge  water  and  thus  are 
returned  to  the  river  in  an  active  form  with  the  return  flow.  These 
dissolved  substances  are  not  subject  to  settling.  The  analytical 
model  can  be  used  to  calculate  the  concentration  of  a  dissolved  sub¬ 
stance  at  any  point  in  the  suspended  solids  plume  with  minor  modifica- 


Equation  3.14  is  the  general  solution  descr-ib’ng  dispersion  and 
settling  in  a  river.  The  exponential  term  describts  the  settling  while 
C'  (y,x)  describes  the  dispersion.  The  equation  that  describes  the 


dispersion  of  a  dissolved  substance  is  then: 
C  (y,x)  =  C‘  (y , x ) 


(5.15) 


LONGITUDINAL  DISTANCE,  X  (meter*) 
4 - DOWNSTREAM 


Figure 


LATERAL  DISTANCE,  Y  (meter*) 


5- ^ 5-  Depth  of  sediment  deposited  from  plume  at 
Rock  Island,  Illinois,  mm. 
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LATERAL  DISTANCE,  Y  (meters) 


Figure  5-16.  Depth  of  sediment  deposited  from  plume  at 
Keithsburg,  Illinois,  mm. 
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The  particular  solution  to  Equation  5.15  is  Equation  3.22.  Examina¬ 
tion  of  Tables  5-3  and  5-4  shows  that  all  of  the  terms  of  Equation  3.22 

Qb 

are  known.  It  becomes  an  easy  matter  to  calculate  *-r—  C'  (y»x)  at 

4oLo 

the  specific  points  (y,x)  in  the  plume.  Multiplying  this  term  by  the 
concentration  of  dissolved  substance  at  the  head  of  the  plume  will 
yield  the  concentration  of  the  dissolved  substance  at  the  several 
points  in  the  plume. 

An  alternate  method  of  arriving  at  the  concentration  of  a  dis¬ 
solved  substance  in  the  plume  is  calculation  of  the  dilution  volume. 

The  dilution  volume  is  defined  as  the  number  of  volumes  of  river  water 
added  to  one  volume  of  water  at  the  head  of  the  plume  to  arrive  at  the 
concentration  in  the  plume.  The  dilution  volume  is  calculated  as: 


Dilution  Volume  = 


C 


o  r 

where  C.  ='vo  o  ,  the  concentration  at  the  head  of  the  plume. 

% 


(5.16) 
The  con¬ 


centration  of  a  dissolved  substance,  C’,  is  then  calculated: 


r 1  =  — 


(5.16a) 


Dilution  Volume 

Tables  5-7  and  5-8  show  the  method  of  calculation  of  the  dilution 
factor  for  the  plumes  at  Rock  Island  and  Keithsburg,  respectively.  It 
can  be  seen  from  Figure  5-17  that  the  concentration  of  a  dissolved 
substance  along  the  shoreline  is  the  same  as  the  concentration  of  that 
substance  at  the  head  of  the  plume,  after  initial  mixing  has  taken 
place  between  the  discharge  water  and  the  river.  This  indicates  that 
settling  is  the  prime  mechanism  operating  to  reduce  the  suspended 
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Table  5-7.  Calculation  of  dilution  volume  for  dissolved  substances,  Rock  Island. 


Dilution 


Table  5-7  (continued). 
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Table  5-7  (continued). 
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Table  5-8  (conti 


solids  concentration  and  that  dilution  played  only  a  small  part  at 
Rock  Island  after  the  initial  mixing.  It  also  indicates  that  a  toxi¬ 
cant  desorbing  from  the  sediment  and  re-entering  the  river  would  be 
in  its  greatest  concentration  along  the  shoreline.  Figure  5-18  for 
Keithsburg  shows  a  much  greater  degree  of  dilution  than  occurs  at 
Rock  Island.  This  is  primarily  because  of  the  much  narrower  plume 
source  with  lower  flow  and  momentum  after  having  traversed  the  disposal 
island. 

The  final  physical  parameter  given  by  the  model  is  the  maximum 
length  of  the  plume.  The  limits  of  the  plume  are  defined  as  the  point 
where  the  suspended  solids  concentration  in  the  plume  is  no  longer 
distinguishable  from  the  ambient  suspended  solids  concentration.  A 
practical  value  for  defining  the  limits  of  a  suspended  solids  plume 
in  the  Mississippi  River  might  be  10  mg/1  above  ambient.  To  find  the 
maximum  length  of  the  plume,  the  model  can  be  solved  for  the  center- 
line  concentration  at  various  values  of  x  until  the  distance  where 
the  suspended  solids  concentration  is  10  mg/1  is  found. 


Ill 


CHAPTER  VI 
WALDEN  PLUME  MODEL 


Walden  Plume  Model 

One  of  the  computational  models  utilized  to  investigate  the 
turbidity  plume  caused  by  dredge  disposal  is  the  Walden  Plume  Model 
presented  by  Wechsler  and  Cogley  (1977).  The  turbidity  plume  model 
was  developed  to  predict  the  suspended  sediment  concentration  down¬ 
stream  from  a  line  source  in  open  water.  The  model  uses  sedimentation 
data  obtained  from  jar  tests  and  hydraulic  data  based  on  simplifying 
assumptions  of  unidirectional  constant  flow,  essentially  infinite 
width,  constant  depth,  and  infinite  length. 

The  mathematical  model  is  a  material  balance  among  the  sediment 
transport  mechanisms  of  (1)  downward  transport  by  settling  with  ulti¬ 
mate  sediment  removal  by  deposition  on  the  bottom,  (2)  upward  transport 
by  vertical  eddy  diffusion  in  the  direction  of  decreasing  concentration 
gradient,  (3)  lateral  dispersion  by  eddy  diffusion,  and  (4)  downstream 
dispersion  by  both  bulk  advection  and  eddy  diffusion. 

The  differential  equation  expressing  the  material  balance  down¬ 
stream  from  a  dredging. site  may  be  expressed  as: 

h  (uc)  +  h  f(w)dw *  ■ 


3_  /r  3_ 

ax  '  x  ax'  ■  ay 


(E  — ) 

y  3y 


■L-  (e  l£)  =  o 

az  [  Z  dZ  u 


(6.1) 
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in  which 

x  =  downstream  distance,  m 

y  =  vertical  distance,  m 

z  =  lateral  distance,  m 

u  =  stream  velocity  at  any  point,  m/sec 

3 

c  =  sediment  concentration,  kg/m 
W  =  settling  velocity,  m/sec 

f (W)  =  settling  velocity  frequency  distribution  (sediment  mass/W) 
vs.  W 

E  ,E  ,E  =  eddy  diffusi vities  in  x,  y,  and  z  directions,  respect- 

As  y  L 

2 

ively,  m  /sec 

In  the  derivation  of  this  equation  it  was  assumed  that  the  flow  is 
steady,  uniform,  and  fully  turbulent,  and  that  eddy  diffusion  can  be 
characterized  by  Fick's  Law  with  eddy  diffusion  coefficients. 

To  apply  the  equation  to  the  plume  model  other  assumptions  must 
be  made  which  are  listed  below. 

(1)  Eddy  diffusion  in  the  downstream  direction  is  negligible  com¬ 
pared  to  the  other  diffusive  transport  terms;  i.e., 

&  (Ex  If)  ■  0  (6-2) 

(2)  The  fully  turbulent  velocity  profile  is  flat,  and  it  can  be 
assumed  that  u  is  constant  and  equal  to  the  mean  velocity, 

U ;  i.e., 

ih<"c>'u§f  (e-3> 

(3)  The  equation  relating  eddy  diffusivity,  E^,  and  vertical 


position  in  the  flow  can  be  derived  by  classical  sediment 
transport  mechanisms  as  described  by  Wechsler  and  Cogley 
(1977). 

Ey  =  0.02  Uy(l  -  y/h)  (6.4) 

in  which  h  =  channel  depth,  m. 

(4)  The  lateral  eddy  diffusivity,  Ez>  is  approximately  constant 
and  can  be  expressed  in  terms  of  the  maximum  value  of  Ey  (at 
mid-depth)  as 

Ez  a  2.2  (.005  Uh)  (6.5) 

Based  on  these  assumptions.  Equation  (6.1)  becomes 


u  If  -  fy  I  /«  f<U>dW  f  -  h  'Ey  fy  >  -  Jz  <EZ  If)  <6-6> 

with  Ey  and  Ez  given  by  Equations  (6-4)  and  (6.5).  The  first  term  of 
Equation  (6.6)  represents  downstream  advection,  the  second  terms 
accounts  for  vertical  sedimentation  (or  settling),  and  the  last  two 
terms  represent  eddy  diffusion  in  the  vertical  and  lateral  directions, 
respectively.  The  integral  in  the  settling  terms  accounts  for  the 
range  of  settling  velocities  of  the  sediment  components  and  must  be 
evaluated  over  the  entire  range  of  settling  velocities.  However,  for 
nonflocculent  sediment,  the  settling  velocity  of  each  particle  is 
invariant  in  time  and  space.  Thus,  the  settling  term  may  be  replaced 
by  the  simpler  form  W  3c/3y,  and  the  equation  which  becomes 


. -  3c  .  i.  3c  3  / c  3c \  3  , r  3C\  _  n 
U  3x  +  W  37  -  37  (Ey  3y}  "  31  {Ez  3?}  ~  0 


(6.7) 


must  be  solved  for  each  sediment  size  present.  The  results  are  then 
superimposed  to  obtain  total  concentration  at  each  point  downstream 
from  the  source. 


If  the  sediment  is  flocculating,  the  problem  becomes  more  diffi¬ 
cult,  and  some  simplifying  assumptions  become  necessary.  This  simp¬ 


lification  is  discussed  in  detail  by  Wechsler  and  Cogley  (1977). 

Equation  (6.7)  is  then  to  be  solved  for  each  settling  -  velocity 
fraction  in  the  sediment.  The  total  suspended  -  sediment  -  concentra¬ 
tion  profile  is  obtained  by  adding  the  concentration  profiles  for  each 
sediment  fraction. 


Numerical  Solution 

Equation  (6.7)  is  solved  numerically  by  a  f ini te-di fference  method. 
The  boundary  condition  at  the  water  surface  specifies  that  there  is  no 
sediment  flux  across  the  surface;  i.e., 

Ey  fy  +  Wc  =  o  (6.8) 


The  bottom  boundary  condition  states  that  all  sediment  reaching  the 
bottom  is  deposited  and  that  there  is  no  reentrainment.  Therefore, 
at  the  bottom 
dC 


E  —  =  0 
y  3y  u 


(6.9) 


The  initial  concentration  at  the  disposal  site,  x  =  0,  is  known; 
i.e.,  at  any  vertical  level,  y,  -b  £  z  £  b,  and  c  =  cQ.  The  calcula¬ 
tion  proceeds  stepwise  downstream.  Over  each  step,  Ax,  an  implicit 
finite-difference  approximation  of  Equation  (6.7)  is  solved  to  compute 
the  concentration  at  the  end  of  the  step  from  the  concentration  at  the 
start.  The  boundary  conditions  at  the  top  and  bottom  surfaces  enter 
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It  is  important  to  note  that  the  longitudinal  diffusion  terms, 
3/3x  (E  3c/3x),  is  expressed  at  x  +  Ax  rather  than  at  x  giving  rise 

A 

to  an  implicit  finite-difference  scheme.  This  approximation  is 
necessary  to  allow  the  use  of  large  values  of  Ax  without  generating 
numerical  instabilities.  The  equation  is  solved  at  N  levels  in  the 
vertical  direction  where  N  =  h/Ay.  The  implicit  system  requires  the 
solution  of  N  coupled  equations  at  each  downstream  step. 

The  effect  of  lateral  diffusion  can  be  found  by  using  the  ana¬ 
lytical  solution  of  the  diffusion  equation  together  with  the  numerical 
solution  previously  described.  The  lateral  diffusion  may  be  described 


1_  (uC)  =  <L  (e  1C) 
3x  VUL;  3z  k  Z  3z' 


(6.10) 


in  which  C  represents  the  concentration  for  the  two-dimensional  problem 
at  each  vertical  level.  The  solution  of  Equation  (6.10)  is: 


4ttxE7  u  r 

c(x,z)  =  nj 


,z  -  v;  u 
4xE_ 


(6.11) 


in  which  v  is  a  dummy  variable  representing  distance  within  the  plume. 
A  transformation  of  variables  relates  this  expression  to  the  error 


function  (erf).  If  Ez  is  taken  as  constant,  and 

■(.z.^_v)2u  „  l£ 

4xEz  2 


Equation  (6.11)  becomes- 

cix'z)  \~k  / 


-Y2/2  j 

e  1  dy 


(6.12) 


(6.13) 


A  transformation  of  the  limits  of  integration  yields  the  error  function 
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for  which  a  solution  is  well  known.  Thus,  an  analytical  solution  can 
be  attained  for  the  lateral  diffusion. 

The  combination  of  the  analytical  solution  for  the  lateral  dif¬ 
fusion  problem  and  the  finite-difference  solution  of  the  sedimenta¬ 
tion  and  vertical  diffusion  problem,  as  described  by  Wechsler  and 
Cogley  (1977),  yields  a  good  approximation  to  the  three-dimensional 
concentration  field,  c(x,y,z). 

Computer  Program 

A  FORTRAN  IV  computer  program  to  predict  the  three-dimensional 
sediment  plume  was  presented  by  Wechsler  and  Cogley  (1977).  Appropri¬ 
ate  revisions  were  made  by  the  present  authors  to  solve  the  problem 
of  interest  herein  using  a  Control  Data  Corporation  (CDC)  CYBER  71 
computer  system.  The  program  is  listed  in  Appendix  E. 

The  use  of  the  program  is  discussed  next.  Input  data  include  the 
stream  velocity,  U,  stream  depth,  H,  sediment  settling  velocity  dis¬ 
tribution  (given  as  the  number  of  sediment  fractions,  NSEDF,  and  the 
concentration,  CO,  and  settling  velocity,  W,  of  each),  the  initial 
discharge  half-width,  XL,  and  two  computational  parameters  -  the 
number  of  downstream  steps,  NSTEP,  and  the  size  of  the  computational 
steps  in  the  lateral  (z)  direction,  DELZ.  The  longitudinal  step  size, 
DELX,  is  taken  as  a  constant;  it  is  defined  in  a  substitution  state¬ 
ment  in  the  program.  The  vertical  step  size,  DELY,  is  computed  from 
the  stream  depth,  H,  and  the  number  of  vertical  steps,  XN. 

Each  sediment  factor  is  analyzed  separately,  and  all  fractions  are 
combined  to  show  the  three-dimensional  sediment  plume.  The  first 


input  value  which  is  entered  is  the  number  of  sediment  parameters 
needed.  The  next  data  which  are  read  in  are  U,  W,  H,  CO,  NSTEP,  XL, 
DELZ.  The  values  change  for  each  sediment  fraction,  and  they  are 
entered  separately  as  the  program  considers  each  fraction. 

After  calculating  certain  constants  to  be  employed  during  the 
program  execution  and  after  computing  the  lateral  eddy  diffusivity 
from  Equation  (6.5),  the  lateral  diffusion  is  found  by  solving  Equa¬ 
tion  (6.13),  the  transformed  version  of  Equation  (6.11).  The  program 
then  performs  a  finite-difference  solution  of  the  longitudinal  dif¬ 
fusion  equation  at  each  level  in  the  vertical  direction. 

The  analytical  and  numerical  results  are  combined  (as  described 
by  Wechsler  and  Cogley,  1977)  by  rewriting  the  diffusion  equation 
symbolically  as 

u|^=(Ly  +  Lz)c  (6.14) 

in  which  Lzrepresents  lateral  diffusion  and  Ly  represents  vertical 
diffusion  and  sedimentation. 

If  C(x,z)  represents  the  analytical  solution,  and  C'(x,y)  repre¬ 
sents  the  finite-difference  solution  of  the  two-dimensional  problem 
ignoring  Lz>  the  required  solution  which  satisfies  Equation  (6.14), 
after  matching  the  initial  conditions,  is 

c(x,y,z)  =  C(x,z)C (x,y)  (6.15) 

Since  Ly  and  Lz  affect  only  C'(x,y)  and  C(x,z),  respectively.it  is 
possible  to  compute  the  analytical  and  numerical  results  separately  and 
combine  them  to  obtain  a  valid  numerical  approximation  of  the  three- 
dimensional  concentration  field. 
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Resul ts 

The  program  output  consists  of:  (1)  the  vertical  distribution  of 
sediment,  in  the  absence  of  lateral  spreading,  downstream  from  the 
source  for  each  sediment  fraction;  (2)  the  summation  of  all  the  verti¬ 
cal  slices  for  all  sediment  fractions;  (3)  the  lateral  spreading 
coefficients;  and  (4)  horizontal  slices  through  the  three-dimensional 
plume  at  preselected  depths. 

Sample  results  are  given  for  the  dredge  disposal  plume  in  the 
Mississippi  River  at  two  sites  -  Rock  Island,  Illinois,  River  Mile  (RM) 

482,  and  Keithsburg,  Illinois,  RM  428. 

Rock  Island,  Illinois 

At  Rock  Island  the  sediment  is  assumed  to  consist  of  25  percent 

silt,  30  percent  clay,  and  45  percent  sand.  The  settling  velocity  of 

-4  -6 

the  silt  is  taken  as  3  x  10  m/sec  and  of  the  clay  as  3  x  10  m/sec 

which  are  in  agreement  with  pipette  measurements  by  Birks  (1980).  The 

settling  velocity  of  the  sand  is  0.012  m/sec.  The  channel  depth  is 

approximately  2  m,  and  the  stream  velocity  is  0.4  m/sec.  The  initial 

width  of  the  sediment  disposal  plume  is  taken  as  25  m,  and  the  initial 

concentration  is  125  mg/1.  The  lateral  dispersion  coefficient,  Ez,  is 

computed  from  Equation  (6.5),  and  it  has  a  magnitude  of  0.0088  m  /sec. 

The  next  three  figures  give  the  vertical  distribution  of  each 
sediment  component,  in  percent  of  total  concentration,  downstream  from 
the  disposal  site.  Figures  6-1,  6-2,  and  6-3  show  (in  tabular  form), 
the  hypothetical  distribution  in  the  absence  of  lateral  spreading  of 
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Figure  6-1.  Vertical  concentration  distribution  of  silt  downstream  from 
source  (no  lateral  spreading)  -  Rock  Island. 
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silt,  clay,  and  sand,  respectively.  It  is  interesting  to  note  that 
the  heavy  sand  particles  are  seen  to  settle  out  of  the  plume  in  a 
relatively  short  distance  downstream  from  the  source.  The  lighter 
silt  and  clay  remain  dispersed  in  the  flow  and  are  transported  a  much 
greater  distance  downstream.  The  concentration  of  the  clay  is  seen  to 
change  very  little  downstream  from  the  source,  since  it  becomes  a 
colloidal  suspension.  The  distribution  of  all  the  sediment  in  the 
plume,  i.e.,  the  summation  of  the  vertical  distributions  of  each  sus¬ 
pended  sediment  fraction  is  shown  in  Figure  6-4. 

Lateral  spreading  factors  are  given  in  Figure  6-5.  These  factors 
are  then  applied  to  the  vertical  concentration  profiles,  and  the 
three-dimensional  plume  is  calculated  as  horizontal  sections  at 
specified  depths.  Figures  6-6  through  6- 1 0  show  the  horizontal  distri¬ 
butions  of  sediment  at  the  surface  and  at  depths  of  0.4  m,  0.8  m, 

1.2  m,  and  1.6  m.  These  figures  give  a  good  picture  of  the  three- 
dimensional  plume  in  tabular  form. 

because  the  Walden  Plume  Model  is  designed  to  simulate  open-water 
disposal,  its  use  in  the  bank-disposal  problem  considered  herein  re¬ 
quires  the  horizontal  concentration  distribution  to  be  folded  about 
the  line  taken  as  the  river  bank.  In  other  words,  a  reflection  prin¬ 
ciple  is  used  to  simulate  the  river  bank.  This  reflection  is  illustrat¬ 
ed  in  Figure  6-11  with. the  horizontal  concentration  distribution  at 
the  depth  of  1.2  m.  The  open-water  plume.  Figure  6-9,  is  folded  about 
a  line  at  the  edge  of  the  line  source,  at  z  =  -12.50  m,  to  approximate 
the  river  bank.  To  determine  the  distribution  for  the  bank  disposal 
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Figure  6-6.  Sediment  concentration  distribution  in  horizontal 
plane  at  the  surface  -  Rock  Island. 
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Figure  6-7.  Sediment  concentration  distribution  in  horizontal 
plane  at  depth  of  0.4  m  -  Rock  Island. 
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Figure  6-8.  Sediment  concentration  distribution  in  horizontal 
plane  at  depth  of  0.8  m  -  Rock  Island. 
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Figure  6-9.  Sediment  concentration  distribution  in  horizontal 
plane  at  depth  of  1.2  m  -  Rock  Island. 
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Figure  6-10.  Sediment  concentration  distribution  in  horizontal 
plane  at  depth  of  1.6  m  -  Rock  Island. 
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Open  Water  Disposal  (see  figure  6-9) 
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Figure  6-11.  Illustration  of  reflection  principle 
for  bank  disposal . 
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(lower portion  of  Figure  6-11),  the  distribution  values  for  the  open- 
water  disposal  (upper  portion  of  Figure  6-11)  along  the  folding  axis 
are  doubled,  and  the  values  below  the  folding  axis  are  determined  by 
adding  the  corresponding  values  of  the  open-water  plume  above  and  be¬ 
low  the  folding  axis. 

A  graphical  presentation  of  the  numerical  simulation  of  the  dis¬ 
posal  plume  at  Rock  Island  is  shown  in  Figure  6-12.  Several  sediment 
concentration  isopleths  are  drawn  at  the  depth  of  1.2  m.  This  depth 
is  of  interest  because  it  is  the  depth  of  the  mean  concentration  of 
the  suspended  solids,  and  it  is  approximately  equal  to  the  depth  where 
field  measurements  were  taken.  Although  field  measurement  data  are 
not  shown  in  this  figure,  the  results  obtained  from  the  Walden  Plume 
Model  are  in  close  agreement  with  the  field  measurements  and  the  re¬ 
sults  given  in  Chapters  IV  and  V.  The  observed  lateral  soread  was 
somewhat  greater  than  that  shown  in  Figure  6-12. 

Another  indication  of  the  plume  orientation  is  given  by  a  plot  of 
the  line  of  maximum  sediment  concentration  at  the  depth  of  1.2  m  with 
distance  downstream  from  the  source.  This  simulation  is  shown  for 
Rock  Island  in  Figure  6-13.  The  field  data  from  Figure  2-8  are  shown 
in  this  figure  as  a  verification  of  the  Walden  Plume  Model.  The  agree¬ 
ment  of  the  model  prediction  and  the  data  is  quite  good.  The  problem 
downstream  from  about  T75-350  m  was  caused  by  the  disposal  operation 
pumping  pure  water  while  the  field  measurements  were  taken  as  explained 
earlier.  Field  observations  indicate  a  larger  lateral  spread  than 
calculated  in  Figure  6-12,  so  the  lateral  dispersion  coefficient  should 
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Figure  6-12.  Numerical  simulation  of  disposal  plume  at  depth  of 
1.2  m  -  Rock  Island,  Illinois. 


maximum  sediment  concentration  at  depth  of  1 
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have  been  greater  than  88  cm  /sec.  It  is  noted  that  the  maximum  con¬ 
centration  levels  out  at  about  150  m  downstream  from  the  source  to 
approximately  66  mg/1.  This  trend  indicates  that  in  the  first  150  m 
downstream  from  the  source,  the  heavy  sand  particles  settle  out  of  the 
plume,  and  the  light  silt  and  clay  particles  remain  suspended  for  a 
long  distance.  A  comparison  with  Figures  6-1,  6-2,  and  6-3  verifies 
this  observation. 

One  of  the  important  variables  in  a  turbidity  study  is  the  set¬ 
tling  velocity  of  the  sediment  particles.  The  settling  velocity  is 
related  to  the  nominal  size  of  the  sediment  particle  which  may  be 
determined  by  a  pipette  and  visual  accumulation  tube  analysis.  Ob¬ 
viously,  the  heavy  particles  have  a  high  settling  velocity  and  settle 
out  of  the  turbidity  plume  first.  The  effect  of  the  sand  settling  vel¬ 
ocity  on  the  maximum  concentration  distribution  at  the  depth  of  1.2  m 
is  shown  in  Figure  6-14  in  which  results  are  given  for  two  different 
settling  velocities  of  sand.  As  expected  the  plume  with  the  heavier 
sand  reaches  its  asymptotic  concentration  first. 

The  effect  of  the  amount  of  sand  in  the  sediment  on  the  turbid¬ 
ity  plume  is  shown  in  Figure  6-15.  The  maximum  solids  concentration  at 
a  depth  of  1.2  m  downstream  from  the  source  is  shown  for  sediments 
with  different  amounts  of  sand.  The  stream  velocity  and  the  sand 
settling  velocity  are  held  constant,  and  the  amount  of  sand  in  the  sed¬ 
iment  is  varied  from  2  percent  to  45  percent.  As  expected,  the  more 
sand  there  is  in  the  sediment,  the  lower  the  suspended  solids  concen¬ 
tration  becomes,  and  the  sooner  the  concentration  levels  out 
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downstream  from  the  source.  After  the  sand  is  settled,  first  the  silt 
and  next  the  clay  fractions  can  be  expected  to  settle,  but  at  distances 
on  the  order  of  ^u/w. 

Figure  6-16  shows  the  effect  of  river  velocity  on  the  suspended 
solids  concentration.  It  is  seen  that  higher  stream  velocities  tend 
to  keep  more  solids  in  suspension,  since  over  a  given  distance  the 
sediment  has  less  time  to  settle  out.  Thus,  a  higher  solids  concen- 
traction  is  maintained  at  any  distance  downstream  from  the  source  by 
a  higher  river  velocity. 

Keithsburq,  Illinois 

The  sediment  at  Keithsburg  is  assumed  to  consist  of  67  percent 
silt,  31  percent  clay,  and  2  percent  sand.  The  settling  velocities 
of  the  sediment  components,  channel  depth,  and  river  velocity  are  the 
same  as  those  for  Rock  Island.  However,  the  initial  width  of  the  dis¬ 
posal  plume  is  3  m,  the  initial  concentration  is  75  mg/1 ,  and  the 

? 

lateral  dispersion  coefficient  is  taken  as  Ez  =  0.03  m  /sec. 

The  two-dimensional  (vertical)  concentration  distributions  in  the 
absence  of  lateral  spreading  are  shown  in  tabular  form  in  Figures  6-17, 
6-18,  and  6-19  for  silt,  clay,  and  sand,  respectively.  The  summation 
of  the  vertical  distributions  of  each  fraction  is  given  in  Figure  6-20, 
and  the  lateral  spreading  factors  are  shown  in  Figure  6-21. 

The  horizontal  distributions  of  sediment  at  the  surface  and  at 
depths  of  0.4  m,  0.8  m,  1.2  m,  and  1.6  m  are  given  in  Figures  6-22 
through  6-26.  The  dredge  disposal  at  Keithsburg  also  was  on  the  river 
bank,  so  the  reflection  principle  must  be  applied  in  the  interpretation 
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Figure  6-23.  Sediment  concentration  distribution  in  horizontal 
plane  at  depth  of  0.4  m  -  Keithsburg. 
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Figure  6-24.  Sediment  concentration  distribution  in  horizontal 
plane  at  depth  of  0.8  m  -  Keithsburg. 
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Figure  6-26.  Sediment  concentration  distribution  in  horizontal 
plane  at  depth  of  1.6  m  -  Keithsburg. 


of  the  horizontal  distributions.  A  graphical  presentation  of  the  hor¬ 
izontal  distribution  at  the  depth  of  1.2  m  is  given  in  Figure  6-27. 

It  is  interesting  to  note  that  due  Lo  the  low  concentration  of  sand 
and  the  relatively  high  dispersion  coefficient,  the  turbidity  plume 
spreads  very  rapidly. 

Figure  6-28  presents  the  variation  of  the  maximum  sediment  concen¬ 
tration  with  distance  downstream  together  with  the  field  data.  The 
agreement  of  the  model  prediction  and  the  field  measurements  is  quite 
good. 

A  complete  set  of  horizontal  sediment  distributions  at  the  depth 
of  1.2  m  is  given  in  tabular  form  in  Appendix  F.  In  these  studies,  the 
stream  velocity,  the  percentages  of  sand,  silt  and  clay  in  the  sediment, 
and  the  fall  velocity  of  the  sand  were  varied.  For  some  of  these 
studies,  the  lateral  dispersion  coefficient  and  the  initial  sediment 
concentration  also  were  varied.  As  mentioned  earlier,  the  reflection 
principle  must  be  applied  to  interpret  these  results  for  bank  disposal. 

The  distributions  are  separated  according  to  sediment  composition. 
For  each  sediment  composition,  sand  fall  velocities  of  0.007  m/sec, 

0.012  m/sec,  and  0.015  m/sec  are  studied.  For  each  fall  velocity, 
river  velocities  of  0.2  m/sec,  0.4  m/sec,  and  0.8  m/sec  are  considered. 
In  some  cases  the  magnitude  of  the  lateral  dispersion  factor  and/or 
the  initial  sediment  concentration  also  were  varied. 

It  should  be  noted  that  even  though  an  implicit  finite-difference 
scheme  was  used  to  calculate  the  vertical  sediment  distribution,  some 
of  the  numerical  results  are  seen  to  be  unstable.  See,  for  example, 


LONGITUDINAL  DISTANCE, 


150 


Figure  6-27.  Numerical  simulation  of  disposal  plume  at  depth  of  1.2  m  ■ 
Keithsburg, ‘Illinois. 
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the  results  for  a  sand  fall  velocity  of  0.015  m/sec  coupled  with  a 
stream  velocity  of  0.2  m/sec.  The  combination  of  the  relatively  high 
settling  velocity  of  the  sand  and  the  low  stream  velocity  yields  very 
high  concentration  gradients  which  lead  to  the  instability.  Another 
example  of  instability  is  shown  on  the  last  page  of  this  appendix  where 
a  high  sand  fall  velocity,  Wg  =  0.03  m/sec,  is  coupled  with  a  sedi¬ 
ment  which  is  almost  all  sand  (95  percent)  and  a  relatively  low  river 
velocity,  U  =  0.4  m/sec.  Results  are  not  reliable  in  these  instances. 

This  appendix  can  be  used  most  effectively  to  compare  results  with 
different  values  of  the  governing  variables.  It  can  be  seen  that 
higher  river  velocities  lead  to  higher  sediment  concentrations  down¬ 
stream  from  the  source.  It  also  is  seen  that  higher  sand  fall  velocities 
lead  to  lower  sediment  concentrations  downstream.  Increasing  the 
magnitude  of  the  lateral  dispersion  factor  is  seen  to  increase  the 
lateral  spread  of  the  plume. 

Appendix  F  is  organized  in  three  parts.  The  sediment  composition 
of  the  first  11  simulations  was  45"'.  sand,  25.'  silt,  and  30?i  clay.  This 
composition  was  characteristic  of  a  medium  grain  sand  (^370  y)  which 
was  pumped  onto  the  beach  and  which  immediately  returns  to  the  main 
channel  of  the  river.  This  was  typical  of  the  beach  nourishment  type 
of  disposal  operation  at  Rock  Island  (Figure  4-6).  At  the  shore  line 
the  sediment  size  distribution  was  95"  sand  (Table  4-1),  but  by  the 
time  it  entered  the  river,  it  was  estimated  at  45"  sand,  25r'  silt,  and 
30w  clay.  This  is  the  maximum  percentage  of  sand  that  one  would  ex¬ 
pect  to  measure  in  the  River  near  the  beginning  of  the  plume.  By 
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measuring  the  dredged  material  size  distribution  and  the  river  velocity, 
one  could  find  a  figure  in  the  Appendix  F  of  similar  characteristics 
and  thereby  estimate  the  extent  and  magnitude  of  the  suspended  solids 

O 

plume.  A  lateral  dispersion  coefficient,  E2,  of  0.044  m  /sec  is  sug¬ 
gested. 

The  second  set  of  plots  in  Appendix  F  is  for  a  sediment  composi¬ 
tion  of  2%  sand,  67%  silt,  and  30%  clay,  and  was  characteristic  of 
the  dredge  disposal  operation  at  Keithsburg  (Figure  4-5).  For  this 
case,  the  discharge  ran  across  Willow  Bar  Island  and  lost  all  but  1.8% 
of  its  sand  (see  Table  4-1).  Once  again  if  the  stream  velocity  is 
known  for  a  similar  case,  the  suspended  solids  plume  could  be  chosen 
from  the  11  at  that  sediment  composition  in  Appendix  F.  Also  available 
are  10  simulations  at  an  intermediate  sediment  composition  of  20%  sand, 
25%  silt,  £ind  55%  clay  which  might  be  representative  of  an  island 
disposal  operation  which  rapidly  returns  to  the  channel. 
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CHAPTER  VII 

SUMMARY  AND  RECOMMENDATIONS 

Field  studies  were  conducted  on  three  dredged  sites  (Hannibal, 
Missouri;  Keithsburg,  Illinois  and  Rock  Island,  Illinois).  Turbidity 
and  suspended  solids  measurements  were  taken  0  -  500  m  downstream  from 
the  discharge  site.  Excess  turbidities  in  the  plume  ranged  from  0-33 
nephelometric  turbidity  units  (NTU)  while  excess  suspended  solids  were 
0  -  125  mg/1.  The  plumes  were  shore-attached  and  near  shore  concen¬ 
trations  (centerlines)  were  measurable  as  far  as  500  m  (at  Rock  Island) 
and  were  less  than  75  m  wide. 

Each  dredging  disposal  operation  was  unique  depending  on  whether 
it  was  a  beach  nourishment  or  island  disposal  type  of  operation.  The 
island  disposal  operation  at  Hannibal  was  entirely  impounded  with  no 
return  water  discharge  whatsoever.  It  is  felt  that  the  "worst  case" 
beach  nourishment  disposal  condition  was  monitored  at  Rock  Island. 

Only  if  the  sediment  were  finer  grained  silt  and  clay  would  a  greater 
suspended  solids  plume  develop. 

Channel  maintenance  dredging  at  the  three  sites  did  not  violate 
Iowa  Water  Quality  Standards  of  25  nephelometric  turbidity  units  (NTU). 
Turbidities  greater  than  25  NTU  were  measured  only  at  the  initial  point 
of  runoff  into  the  Mississippi  River  for  the  beach  nourishment  type  of 
dredge  materials  disposal  at  Rock  Island.  Such  short  term  concentrations 
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would  be  within  an  allowable  mixing  zone  of  most  State  Water  Quality 
Standards . 

Three  mathematical  models  were  utilized  to  describe  the  collected 
field  data:  the  Schubel-Carter  (1978)  model,  the  Wechsler-Cogley  (1977) 
Walden  Plume  model,  and  an  analytical  solution  developed  herein.  Pre¬ 
liminary  results  show  that  the  Schubel-Carter  (1978)  nomogram  solution 
is  cumbersome  to  use  for  riverine  conditions  and  involves  a  very  time 
consuming  trial  and  error  technique  to  calculate  the  correct  initial 
suspended  solids  concentration  at  the  point  of  discharge.  The 
Wechsler-Cogley  (1977)  Walden  Plume  computer  model  has  proven  to  have 
several  advantages  over  the  Schubel  and  Carter  approach.  First,  it  is 
possible  to  use  a  plane  source  discharge  which  is  more  realistic  than 
a  line  source  as  in  Schubel  and  Carter  (1978).  Secondly  it  can  handle 
several  size  fractions  easily  and  the  computations  are  quickly  facili¬ 
tated  by  digital  computer.  The  analytical  solution  developed  herein 
utilizes  probability  density  function  tables  and  is  easier  to  understand 
than  the  numerical  solution  of  Wechsler  and  Cogley,  but  it  does  require 
extensive  hand  calculations. 

The  Walden  Plume  model  and  the  analytical  solution  developed  herein 
were  successfully  used  to  simulate  the  shore-attached  centerline  of  the 
dredge  disposal  operations  at  Keithsburg  and  Rock  Island.  It  is 
recommended  that  these?  models  be  used  in  future  modeling  efforts. 

Lateral  concentration  variations  were  not  well  described  due  to 
insufficient  field  data  as  well  as  a  lack  of  knowledge  of  the  lateral 
dispersion  coefficients  under  these  conditions.  It  is  therefore 
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recommended  that  further  studies  be  undertaken  to  better  delineate  the 
lateral  dispersion  phenomena  as  well  as  the  initial  mixing  and  density- 
dependent  settling  at  the  head  of  the  plume.  Furthermore  a  worst  case 
of  beach  (nourishment  disposal  at  a  site  with  silt  or  clay  sediment 
should  be  monitored  if  such  a  situation  arises. 

Each  of  the  two  models  employed  have  relative  advantages.  The 
analytical  solution  was  conveniently  utilized  to  provide  estimates  of 
the  in- situ  dilution  factors  for  dissolved  constituents  as  well  as  the 
expected  rate  and  depth  of  sedimented  material  in  the  River.  The 
Walden  Plume  model  was  used  to  generate  a  range  of  solutions  for  dredge 
disposal  operations  provided  in  Appendix  F.  If  a  planner  or  engineer 
knows  the  grain  size  distribution  of  the  material  to  be  dredged,  the 
approximate  river  velocity,  and  the  mean  depth  of  the  discharge  area,  it 
is  possible  to  locate  a  graph  in  Appendix  F  of  similar  conditions  and 
to  predict  the  extent  and  concentration  of  the  suspended  solids  plume. 
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ANALYTICAL  MODEL  PROGRAM  FOR  HEWLETT-PACKARD  29C 
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TABLE  47 


Mississippi 

10-24  July 

Dilution  Factors 

River  Pre-Dredging 

1979 

for  Mississippi  River-3 

Samples 

River  Miles323. 

.5-4  89  .4 

Pollutant 

Site  Water 

Elutriate 

W  Q 

Di lution 

Concentration 

Concentration 

Criteria 

Factor 

Ca  ^ 

c  y  g/i 

Ds  yql/ 

Vol/Vol 

COD 

27100 

55900 

- 

- 

Unionized 

Nh3~N 

1.35 

40.45 

16.5 

1.58 

PO  -P 

4 

225 

213 

16.5 

1.58 

Arsenic 

2.6 

3.3 

100 

0 

Beryllium 

- 

.  - 

11 

- 

Cadmium 

- 

- 

1.2 

- 

Chromium 

- 

- 

100 

- 

Cppper 

1.6 

2.2 

800 

0 

Iron 

10.7 

154.7 

1000 

0 

Lead 

- 

- 

315 

- 

Manganese 

0.61 

402  .9 

100 

Mercury 

- 

- 

0.1 

- 

Nickel 

- 

- 

100 

- 

Selenium 

- 

- 

250 

- 

Zinc 

8900 

19000 

136 

N.P.* 

•Not  Possible  to  achieve  water  quality  criteria  by  dilution. 
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Manganese  desorbs  occasionally  at  the  sites  sampled  by  the  GREAT  XI 
Simulation  Study  Group  and  in  2  out  of  3  occasions  in  the  U.S.G.S. 
samples.  In  fact,  in  one  instance  it  requires  a  dilution  factor  of 
3.04.  Although  most  of  the  other  pollutants  analysed  desorb  at  the 
two  "dirty"  sites  -  Aspelmeier  Ditch  and  Grey's  Chute,  their  concentra¬ 
tions  hardly  approach  the  Water  Quality  Criteria  Concentrations. 


Dilution  Factor  for  Dissolved  Oxygen 

It  has  been  shown  that  sediment  dredged  from  a  river  bed,  when 
mixed  with  site  water,  can  exert  an  oxygen  demand.  In  fact,  the  oxygen 
demand  could  even  turn  the  dredge  discharge  anoxic,  under  severe 
conditions.  The  following  calculations  were  made  to  determine  whether 
the  dredge  discharge  would  have  been  anoxic  at  the  eight  sites  sampled 
by  the  GREAT  II  Simulation  Study  Group  and  if  so,  the  dilution  volume 
required. 

« 

Taking  the  case  of  Keithsburg  (Musgrove,  1980), 

Total  volume  of  water  and  sediment  pumped 

=  45  ft^/sec  x  3600  sec/hour  x  18.5  hours  x  .02832  m^/ft^ 

=  84875  m^. 

Total  sediment  removed  =  8885  m^. 

Patio  of  water  pumped  to  sediment  removed 


84375  -  3885 
6885 


8.55. 


»  -iss  balance  on  the  oxygen  concentration, 

<  D.C.i ,  ♦  D ( D . 0 . )  =  (1  +  D)  x  C 

*  1  s 

1  */f»d  jxygen  concentration  at  the  exit  of  the 


(28) 


pipe,  mg/1 
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(D.O.)^  =  dissolved  oxygen  concentration  of  the  ambient  site  water, 
water,  mg/1 

C  =  Water  Quality  Standard,  mg/1 

s 

D  =  dilution  factor,  vol/vol 

Assuming  (D.O.)^  =  8  mg/1  and  Cg  =  4  mg/1, 

4  -  (D.O. ) _ 


8-4 
4  -  (D.O.) 


(29) 


To  calculate  the  D.O.  concentration  of  the  dredge  discharge  the 
oxygen  demand  of  the  discharge  is  determined. 

I 

Taking  the  data  for  Grey's  Chute  (Geadelmann,  1979) , 

_  _  ,  ,  506  gm  of  0„  , 

Oxygen  uptake  of  sediment  =  _ 2  m  one  hour. 

3 

m 

(Total  Oxygen  Demand)  60  min. 


Ultimate  Total  Oxygen  Demand 

K  is  the  rate  constant,  K  =  0.16/min.,  t  =  60  min. 


,,  ~Kt, 
(1  -  e  ) 


Ultimate  Total  Oxygen  Demand 


504 


1  -  e 


-0.16  x  60 


=  504  gm/m~ 


Detention  time  in  a  pipe  1000  feet  long  and  20"  in  diameter  with  a 
discharge  rate  of  45  ft3/second. 


A  x  1  x  ( 20/12)  x  1000  „  _  . 

t  *  — - —  =  - — - — -  =  0.808  mm 

Q  45  x  60 

„ _  .  .  ,.  -0.16  x  . 808  .  3 

Oxygen  demand  =  504  (1  -  e  )  =  61  gm/m  . 


Available  oxygen 


8.55  m3  of  water  8  mg  1  g/m3  _  ^  ^2 

-*  x  j  x  nc/1  - 


of  sediment 


m  sediment 
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Since  the  available  oxygen  68.4  gm/m'  is  greater  than  the  oxygen 
uptake  of  the  sediment  61  gm/m^ ,  the  discharge  is  not  anoxic. 

D.O.  concentration  of  discharge  =  88 •  1-  '  —  =  0.85  mg/1. 

o .  j 

Table  48  gives  the  oxygen  demand  at  the  other  sites  using  appro¬ 
priate  K  values  as  calcuated  by  Geadelmann  (Geadelmann ,  1979) . 

However,  even  if  the  discharge  is  deoxygenated ,  only  if  the 
discharge  were  submerged  would  reaeration  be  impossible,  necessitating 
the  calculation  of  dilution  factor.  Table  48.  gives  dilution  factors 
if  the  dredge  discharge  is  submerged. 

D  =  - — 4°'8'5-  =  0.79  for  above  example. 

Kinetic  Model  of  the  Standard  Elutriate  Test 

The  objective  of  the  kinetic  model  is  to  describe  the  dissolved 
concentration  of  the  pollutant  that  has  been  adsorbed  or  desorbed  from 
the  sediment  during  the  course  of  elutriation,  according  to  the  Standard 
Elutriate  Test.  Since  the  Standard  Elutriate  Test  is  a  simulator  of 
dredge  disposal,  the  kinetic  model  would  also  simulate  the  effects  of 
mixing  bottom  sediment  and  site  water  in  the  dredge  line  upto  the  point 
of  discharge  into  the  river  during  dredging. 

Sorption  is  considered  to  be  the  physical  phenomenon  that  predom¬ 
inates  and  the  model  takes  into  account  only  sorption  kinetics. 

Complete  mixing  is  assumed  within  the  control  volume  -  which  is  the 
volume  of  the  beaker  in  which  the  elutriate  test  is  conducted.  The 


driving  force  in  the'kinetics  is  the  difference  between  the  equilibrium 
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TABLE  48 

Dilution  Factors  for  Dissolved  Oxygen 


Site 

TOD 

60 

g  of  02 

3 

m 

TODu 

g  of  02 

3 

m 

TOD.808 

g  of  02 

3 

m 

Available  02 

g  of  02 

3 

m 

Dil. 

Factor 

Vol 

Vol 

(D.0.)2 

22. 

l 

Maquoketa 

272 

273 

21.2 

68.4 

0 

5.5 

Aspelmeier 

392 

392 

47.5 

68.4 

0.39 

2.44 

Montpelier 

96 

97 

7.5 

68.4 

0 

7.12 

Muscatine 

62 

63 

4  .9 

68.4 

0 

7.43 

Keithsburg 

37 

38 

2.9 

68.4 

0 

7.66 

Grey  Chute 

504 

504 

61.1 

68.4 

0.79 

0.85 

Keokuk 

149 

149 

11.6 

,  68.4 

0 

6.60 

Quincy 

56 

57 

4.4 

68.4 

0 

7.48 

TODgo  *  Total  Oxygen  Demand  for  60  minutes. 


TODy  »  Ultimate  Total  Oxygen  Demand. 

1100  ana  m  Total  Oxygen  Demand  for  0.808  minutes. 

•  quo 

(D.O.)2  *  Dissolved  Oxygen  Concentration  of  the  discharge. 


the  control  volume,  then 


CP  *  rn. 

0  0 

However,  mass  must  be  conserved  and  the  total  pollutant  concentra¬ 
tion  in  the  beaker  must  remain  constant.  The  ratio  of  dissolved  and 
adsorbed  pollutant  concentration  will  now  be  dynamic  and  change  to 
approach  an  equilibrium  state  at  which  P  =  r/C. 

Taking  a  mass  balance,  the  total  pollutant  concentration. 


C  «  C  +  r. M  =  C  +  rM 
TOO 
» 

Rearranging  equation  (32) , 

Crt  +  r  M  -  C 
0  0 


M 


(32) 


(33) 


Substituting  for  r  in  (31) 


c  +  r  M  -  C 

§“  “  (JL“ F - PC) 


jp 

i.e.,  —  «  CK  (1  +  PM)  =  K  (CQ  +  rQM) 


(34) 


Comparing  (34)  to 


^  +  P  (x)  y  *  Q  (x) 


the  solution  to  whic'i  is 

fc.dx  ,  $P.dx 
y.e  -JQ-e  +  C 

Hence  solution  to  (34)  is 

Jk(1  +  PM).dT  f  jk(l  +  PM).dT 

C.e  -  \  K(CQ  +  rQM) .e  +  C  (35) 

Using  the  boundary  condition  that,  at  T  =  C,  C  =  CQ  and  simplifying. 
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where  =  dissolved  concentration  of  pollutant  at  time  T  *  0,  jjg/1 


K  =*  sorption  rate  constant,  1/hour 
P  =  partition  coefficient  of  pollutant  at  the  site. 


»*9Ag 


M  =  suspended  solids  concentration,  kg/1 
T  =  time ,  hours 

rQ  =  adsorbed  pollutant  concentration  at  time  t  =  0,  pg/kg 


Experimental  Methods  and  Results 

Two  sites  near  Keokuk  (near  M.P.  355)  -  one  a  main  channel  and  one 
a  slough  (Grey  Chute)  were  sampled  in  November,  1979.  Both  were 
proposed  dredge  sites.  A  Ponar  Dredge  Sampler  was  used  to  collect  all 
the  sediment  samples .  The  Ponar  Dredge  Sampler  is  capable  of  sampling 
to  a  maximum  depth  of  five  inches  and  has  a  volume  of  nearly  three 
liters.  The  site  water  was  collected  about  a  foot  below  the  water 
level  in  5  gallon  Nalgene  containers.  The  sediment  container  and  the 
Nalgene  container  had  been  prewashed  with  acid  and  rinsed  with  deionized 
water. 

The  kinetic  experiments  were  conducted  within  2  days  of  sampling, 
in  an  open  5  gallon  Nalgene  container  using  the  sediment  and  unfiltered 
dredge  site  water  (see  Figure  45)  .  Three  liters  of  sediment  and  12 
liters  of  site  water  were  used  for  the  test  to  give  a  sediment  to 
water  volumetric  ratio  of  1:4.  A  non-metallic  propeller  of  diameter 
20  cm  with  a  stem  of  length  40  cm  powered  by  a  1/80  HP,  5000  PPM  motor, 
with  a  speed  reducer  capable  of  reducing  the  PPM  to  a  fourth  of  the 


normal  speed.  Two  large  course  bubble  stone  diffuser  tubes  wefe  placed 


at  the  sides  of  the  container,  at  the  bottom,  to  aid  in  nixing  and  to 


170 


avoid  anoxic  conditions.  Compressed  air  was  passed  through  a  deionized 
water  trap  before  entering  the  diffuser  tubes  and  slurry. 

o 

Two  tests  were  conducted  at  7.5  C  and  two  at  room  temperature. 

The  motor  was  stopped  during  sampling.  Samples  were  drawn  off  initially 
at  15  minute  intervals  and  then  at  the  second,  fourth,  fifth,  sixth  and 
twenty-fourth  hours.  The  sampling  does  cause  a  reduction  in  volume  of 
the  overlying  water  since  the  samples  taken  do  not  contain  much 
suspended  solids.  The  final  sediment  to  water  ratio  is  approximately 
3.7:1. 

The  samples  were  first  centrifuged  to  reduce  suspended  solids  and 
to  expedite  subsequent  filtering.  A  47  mm  diameter,  5  micron  filter 
was  used  initially  followed  by  a  0.45  micron  filter.  The  filtering 
apparatus  was  entirely  made  of  glass  (millipore)  and  a  vacuum  pump  was 
used  to  facilitate  filtering.  The  samples  were  then  stored  at  4°  C 
until  analysis  for  COD  and  :rH ^  within  a  week. 

Chemical  Oxygen  Demand  was  determined  in  accordance  with  Standard 
Methods  (APHA,  Standard  Methods,  1975).  Ammonia  was  determined 
according  to  procedures  outlined  in  Standard  Methods  (APHA,  Standard 
Methods ,  1975) . 

The  sorption  rate  constant  for  the  four  cases  studied  was 
calculated  by  trial  and  error.  Constant  values  ranging  from  0.01  to 
1 . 3/hour  were  attempted  until  a  rate  constant  of  0.1/hcur  was  found  to 
have  a  satisfactory  fit  to  all  me  experimental  values.  This  value  of 
sorption  rate  constant  was  used  m  mciel  the  experimental  data. 

The  Kinetic  model  was  solved  uni  no  a  FORTRAN  computer  program  on 
tne  CDC  CISC?.  71.  A:  .  nr .  .  :  '■  I.  •  ,,  h.iv«-  L<  <_*ii 
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achieved  after  24  hours,  the  total  pollutant  concentration  was  based  on 
the  24  hour  pollutant  concentration.  The  mass  of  suspended  solids  in 
the  control  volume  was  found  to  be  0.43  kg/1  in  the  case  of  Grey  Chute. 
This  value  of  M  was  used  in  all  the  cases.  The  values  of  partition 
coefficients  were  taken  from  data  already  available  for  the  two  sites 
(Table  35  ) . 

Figure  46  shows  the  model  results  and  the  experimental  results  of 
the  kinetic  experiment  for  COD  at  Grey  Chute  at  4.5°  C.  As  can  be  seen, 
the  model  closely  approximates  the  experimental  values  obtained, 
particularly  until  the  second  hour.  It  should  be  noted  that,  according 
•to  the  model  about  81%  of  the  concentration  change  is  complete  within 
the  first  30  minutes  and  that,  from  the  experimental  values,  95%  of  the 
concentration  change  is  complete  within  the  first  30  minutes. 

Figure  47  gives  the  model  result  and  the  experimental  values  of  the 
kinetic  experiment  for  ammonia  at  Grey  Chute  at  7.5°  c.  Again,  the 
model  closely  approximates  the  experimental  values  throughout  the  test 
period  of  24  hours.  Figure  48  shows  the  model  results  and  experimental 
values  for  COD  at  Grey  Chute  at  room  temperature.  After  the  15th 
minute  the  highly  variable  experiments  do  not  conform  to  the  model 
predictions.  Figure  49  gives  the  model  results  and  experimental  values 
for  OOD  at  Keokuk  at  roan  temperature.  In  this  case  the  experimental 
values  are  quite  close  to  the  model  predictions.  In  the  4  instances 
cited  above  the  experimental  values  showed  that,  on  an  average,  72%  of 
the  concentration  change  is  over  within  the  first  30  minutes.  According 
to  the  model,  on  an  average,  68%  of  the  concentration  change  is  complete 


by  30  minutes . 


K  =  0.1/hr 
Model 

Experimental  Values 


Figure  49.  COD  at  Keokuk  -  Room  Temperature  M.P.  356 


AD” A 086  507  IOWA  INST  OF  HYDRAULIC  RESEARCH  IOWA  CITY  F/«  8/8 

SUSPENDED  SEDIMENT  MODEL I N6  OF  DRED6E-01SPOSAL  EFFLUENT  IN  THE  <— ETC<U> 
MAR  80  J  L  SCHNOOR •  A  R  6IAOUINTA  DACW25-78-C-0073 

UNCLASSIFIED  IIHR-72  NL 
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In  the  case  of  COD  at  Keokuk,  the  ice  bath  malfunction  led  to  a 
temperature  rise  from  8°  to  18°.  A  scrutiny  of  the  experimental  results 
suggested  either  volatilization  of  biological  degradation  of  the  COD. 
Assuming  an  additional  first  order  reaction,  a  degradation  or  volatili¬ 
zation  term  was  added  to  the  original  equation  (equation  31) .  The 
equation  governing  the  kinetics  are  now  as  follows: 

||  =  KM  (r  -  PC)  -  KaC  (37) 

||  *  -K  (r  -  PC)  (38) 

Since  there  is  no  conservation  of  mass,  due  to  either  volatiliza¬ 
tion  or  biological  degradation  of  COD,  the  assumption  of  a  constant 
total  pollutant  concentration  cannot  be  made.  Also,  it  is  not  possible- 
to  calculate  the  total  pollutant  concentration  from  the  24  hour  pollu¬ 
tant  concentration  as  steady  state  is  not  achieved.  The  adsorbed 
pollutant  concentration  used  was  achieved  by  trial  and  error.  The 
degradation  rate  constant  used  was  calculated  through  a  sensitivity 
analysis.  The  simultaneous  equations  (equations  37  and  38)  were  solved 
using  GASP  IV,  which  is  a  Combined  Continuous/Discrete  FORTRAN  based 
Simulation  Language.  The  model  results  and  experimental  values  are 
given  in  Figure  5o  ■ 

The  Kinetic  Model  attempts  to  explain  the  mechanics  of  sorption 
during  the  process  of  elutriation.  The  experimental  results  showed  that 
the  elutriate  test  with  30  minutes  of  agitation  is  a  fair  approximation 
(72%)  of  equilibrium.  The  Kunetic  Model  shows  that  besides  the  sorption 
constant  there  are  two  more  factors  -  the  partition  coefficient,  P,  and 
the  suspended  solids  concentration ,  M,  that  affect  the  dynamics  of 
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elutriation.  The  Kinetic  Model  could  be  used  to  model  the  concentra¬ 
tion  plume  of  a  pollutant  during  disposal,  as  is  shown  subsequently. 


Proposal  of  Plume  Model 

The  model  characterizes  the  pollutant  plume  for  open  water  dredge 
disposal  from  the  pipe  line  of  a  hydraulic  dredge.  The  pollutant  is 
thoroughly  mixed  in  the  pipe  line  before  disposal.  On  disposal,  the 
assumption  made  is  that  there  is  no  dispersion  so  that  the  plume 
concentration  represents  the  worst  case.  The  model  is  one  dimensional 
along  the  center  line  of  the  plume.  The  kinetics  of  sorption  and 
sedimentation  are  assumed  to  predominate  to  the  exclusion  of  all  others. 
Formulation 

The  driving  force  in  the  sorption  kinetics  is  assumed  to  be  the 
difference  between  the  actual  adsorbed  concentration,  r,  and  the 
equilibrium  adsorbed  concentration,  r  .  At  equilibrium, 

r  =  PC  (38) 

eq 

The  model  itself  is  divided  into  two  stages.  The  first  stage  of  model 
describes  the  pollutant  concentration  within  the  pump  and  pipe,  which 
is  similar  to  the  kinetic  model.  The  second  stage  describes  the  pollu¬ 
tant  concentration  and  suspended  solids  concentration  of  the  plume  on 
discharge  from  the  pipe.  The  second  stage  has  a  term  to  take  into 
account  the  settling  of  the  suspended  solids. 

Stage  I: 


dC 

dT 


KM  (r  -  PC) 


(39) 


dr 

dT 


K  (r  -  PC) 


(40) 


SUSPENDED  SOLIDS,  mg/£ 


•  isr*  — 
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Figure  D-2.  Suspended  solids  vs.  turbidity  relationship  for 


Rock  Island  field  data.  All  values  are  expressed  as 
mg/8,  or  NTU  above  ambient. 


APPENDIX  E 

WALDEN  PLUME  MODEL 


COMPUTER  PROGRAM 
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00100  P8DGRAB  PLUHB  (INPUT,  DOTPUT,OOTPLB,TAPE5=r  HPOT,P  APE6=0UTP0r, 

00110*  TAPE8*OOTPLB) 

00120C 

00130; 

00140C  PRO  GRAB  SOURCE  A  LABORATORY  STUDY  CP  THE  TURBIDITY 

001 50C  GENERATION  POP FMTIAL  OP  SEDIBEMTS  TO  EE  DREDGED,  BY 

00160C  B.A.  8RCHSLSR  C  D.R.  C03LEY,  TECH.  REPORT  D-77-14, 

00170C  U.S.  ARBY  ENGINEEER  WATERWAYS  EX  PER  IB  EN  T  STATION,  NOT.  1977 

00180C 

00190C 

00200C  WALDEN  PLUNE  BOOEL 

00210;  ♦*•**#**«*•*.****••. ******************************************* 

00220C  ASSUBPTION  S  INCLUDE  STEADY  STATE,  NO  BCBENTUB  EPPECTS 

00230C  AND  NO  RES  US  PENS IOR  OP  RATERIAL 

00235C 

002 40C  •*•*♦«******•*******••****•*»*********»•*•**•«*•*.*••******».•* 

0024 1C  THIS  PROGRAB  CONTAINS  A  SUBROUTINE  PROP  THE  IHSL  LIBRARY,  A 

00242C  PROPRIETARY  PACKAGE  FRO!*  THE  INTERNATIONAL  B  A  THEN  A  TIC  A  L  G 

00243C  STATISTICAL  LIBRARIES,  INC.,  HOUSTON,  TEXAS.  THIS  ROUTINE  BAY 

00244C  NOT  BE  REDISTRIBUTED  OR  REHOYED  PROB  THIS  SOFTWARE  FDR  OSE  IN 

00 2 4 SC  OTHER  SOFTWARE  D E YE LOPBE NT.  THE  IHSL  ROUTINE  INCLUDED  IS  ERF. 

00246C  *****•#**.«*,**••***,****•***.********.•**•***•..*.*******«*••. 

00247C 

002S0  REAL  IA,BA,HL,IIA 

00260  DOUBLE  PRECISION  B8,HP 

00270  DIHSNSION  C (20 , 5 1 ) , D ( 20 , 20) , I A ( 20, 20) , B A ( 20, 20) , R L ( 400) 

00  280  DIMENSION  R  A  (2  0 , 20)  ,  C  J  (2  0)  .  HL  (4  0  0 )  ,  3  (20  )  ,  BP  (40  3  )  ,  1 1 A  ( 20 , 20) 

00290  DIMENSION  ADELX  (51)  ,  ADELY  (20)  ,‘4  (  21, 51)  ,  AYG  (51 )  ,CSUH  (20, 51 ) 

00300  DIMENSION  A YGSOB (51 ) ,09 (20)  ,IOUT (51 ) 

00310  DATA  CSUB/1020*0./ 

00320  DATA  IA/400*0./,3/400*0./,IIA/400*0./,Z/1050*0./,AYG5UE/51*0./ 

00330C  STATEMENT  FUNCIION  TO  CALCULATE  EDDY  DIFFUSIVITY  AT  ANY  DEPTH 

00340  E  (T)  =0. 02*  D*Y*  (1  .-Y/H) 

00350  DO  10  1=1,20 

00360  IA (I,I)=1. 

00370  XIA (1,1)  =- 1 . 

00300  IP  (I.GT.1)  IIA(I,I-1) = 1. 

00390  10  CONTINUE 

00400  NS2D= 1 

0041 OC  NO.  OF  SSD .  FRACTIONS . 

00420  READ*,NSEDF 

00430C  *********••****•***»»♦**«•**••**••**•*•** **••***••*****♦**• ****** 

0044 OC  U=STRSAB  Y SLOC ITY , B/S EC  W  =  S ETTL IN G  YELOCITY , B/SSC 

00450C  H*STREAM  DEPTH  ,  H  CO=CONCENT RATION  OP  SED.  FRACTION 

00460C  NSTEP=  NO. OF  DOWNSTREAM  STEPS  XL=INITlAL  DISCHARGE  HALP- WIDTH, H 

00470C  DSL  Z=LA  TER  A  L  STEP  SIZS,H. 

004 eo;  *****••*****•*«*••*•***•«**♦•******•**•***••***•***•*••••*******. 

00490C  PROGRAB  CONTROL  IS  TRANSFERRED  HERS  FCR  EACH  SED. FRACTION 

00500  11  IF  (NSBD  .LE.  NS  EOF)  RE A D* , U , W, H, CO, N STEP , XL ,DELZ 

00510  IP ( N SED  .GT.  N SEDP) GO  TO  99 

00520C  XN= NUMBER  OP  DEPTHS 

00530  XN  =  20 

00535  THETA* 1.0 

00540  IP  (H  STEP . GT. 50)  WST8P*50 

00550  NSTEP1= NSPEP*1 

00560  WRITE(8, 1001)  U,W,R , CO, NSTEP, XL, DELZ 

00570  1001  FORMAT (//2 X , *1 NPUP-  *, *U  ■ *, F6. 2, 2X, * W  =*,FB.6,2X, 

00580*  *H  =*,F5. 1, 2X,*CO  =• ,F 6 . 2 ,2 X, * NSTEP  =*,I5,2X, 

00590*  *XL  = • , F5 . 2/, *  DELZ  =*,F6.2//) 

00600  N  =  XN 
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0061 OC  LONGITUDINAL  STEP  SIZE.,.,...,... 

00620  DELX-10. 

00630C  VERTICAL  STEP  SIZE . , . 

00640  DEL Y=H/XN 

00650  CONZ-2.2 

00660  DELI2=1./(DELT-DELT) 

00670  A YG SUB  (1 )  =  AVGSU H  ( 1)  -CO 

00680  INC= 1 

00690  DO  501  J=1 , NSTEP1 

00700  ADELX(J) =(J-1) -DELX 

00710  AVG(J)=0. 

00720  501  CON  TIN OE 

00730  DO  500  1-1, N 

00740  ADELT(I) =  (I-1) -DELT 

00750  IP  (  (1-1 )  *  DELZ.  LE.  X  L)  Z(I,1}=1. 

00760  500  CONTINUE 

00770C  CALCULATE  LATERAL  EDDY  DIFFOSIVITY, EQ. (5) 

00780  EZ=0.005*H«0*COSZ 

00810  FOREZ-4 . *EZ 

00820C  CALC0LAT2  CONCENTRATION  DUE  TO  LATERAL  DISPERSION 

00830  DO  690  J=  1 , NST EP 

00840  F0BEX=SQBT{ADELX(J-1) *POREZ) 

00850  DO  630  IZ-1,21 

00860  AZ=  (IZ- 1)  *DELZ 

00870  TOP= (AZ-XL) /FOREX 

00980  ET=EB  F ( TOP) 

00990  BOT=  (AZ-XL) /FOREX 

00900  EB  =  £R  F (BOT) 

00910  Z(IZ,J-1)*0.S- (BT-EB) 

00920  680  CONTINUE 

00930E  START  FINITE  DIFFERENCE  SOLUTION . 

00940  DO  100  1=1,20 

00950  C(I,1)=C0 

00960  100  CONTINUE 

00970  A V3  (1 ) =C0 

00980  0(  1,  1)  —B(1. 5-DELI) 

00990  D  (1 ,2) =E  (1. 5-DELY) 

01000  N 1 =N-  1 

01010C  CALCULATE  EDDY  DIFFUSIVITY  AT  VARIOUS  DEPTHS 

01020  DO  200  1=2, Ml 

01030  XI- (2-1-1) -. 5-DELY 

01040  X2= ( 2*1 ♦ 1) -.5-DELY 

01050  D  (I ,I-1)  =  B  (XI) 

01060  D(I,I)  — E(X1)-E(X2) 

01070  D  (I  ,1-1 ) -S  (X2) 

01080  200  CONTINUE 

01090  D  (N,N)*-E ( (XN-.5) -DELI) 

01100  D(N  ,N-  1)  =-D (N,N) 

OHIO  DO  300  1=1, N 

01120  DO  300  J=1,N 

01130  D  (I  ,  J)  =  DSLI2-D  (I,  J) 

01140C  CALCULATE  COEFFICIENTS  OF  SISTER  OF  EC.  RESULTING  FROM 

01150C  FINITE  DIFFERENCE  SOLUTION 

01160  NA( I,J) -U/UELX-IA (I , J) -THETA-D (I ,J) 

01170  BA (I , J) -O/DELX -I A ( I , J)  ♦ { 1. -THETA) -D  (I, J ) -W-IIA (I, J) /DELI 

01180  300  CONTINUE 

01190  CALL  ARRAY  (2, N,N, 20,20, RL.RA) 

01200  CALL  ARBAT (2,N , N ,20,20,RL,BA) 

01210  DO  400  J*  1  , NST  EP 

01220  DO  405  1*1, N 

01230  CJ(I)-C(I,J) 
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01200  CS0H(X,J)  =csm  {I  ,J)  ♦2J{I) 

01250  005  CONTINUE 

01260C  FINAL  SOLUTION  BY  GOBBI NATION  OF  ANALYTIC  AND  NOB.  SOLUTIONS 

01270  CALL  GBPRD  (RL,CJ,B, N,N, 1) 

01280  DO  006  1=1,000 

01290  BP  (I)  =HL  (I) 

01300  406  CONTINUE 

01310  DO  407  1=1,20 

01320  B8(I)=B(I) 

01330  407  CONTINUE 

O1340C  SOLVE  SYSTEB  OP  EQUATIONS . 

01350  CALL  DGELG  (B9, BP, N,  1, .000000001, IER) 

01360  IF  (IER.GT.O)  WHITE (6,4747)  IER.J 

01370  4747  FOR  BAT (*  LOSS  OF  SIGNIFICANCE  AT  PIVOT  *,I3,*  IN  STEP  *,I3) 
01380  DO  399  K=1,N 

01390  C  (K,  J’  1 )  =B8  (!C) 

01400  AVG(J*1)  =B8(K)  »AVG(JM) 

01410  399  CONTINUE 

01420  AVG(J  +  1)  =A  VG  (J’l)  /XN 

01430  AVGSUB  (J*1)=AVGSUB(J»1)  ’AVG(J*1) 

01440  400  CONTINUE 

01450  WRITE (8 ,5001)  ( A  DELL  ( I) , 1= 1 , NSTEP 1 , 5) 

01460  5001  FORHAT( 1 1X,20F6. 1 ,/) 

01470  DO  410  1=1 ,N 

01480  CSUB(I,NSTEP1) =GSU!i (I,NSTEP1)*C (I,SSTEP1) 

01490  WRITE  (8,5000)  ADELY  ( I)  ,  (C  (I,  3)  ,  J  =  1,  N  STEP  1 , 5) 

01500  5000  FOHBAT(3X,F5.2,2X,20F6.2) 

01510  410  CONTINUE 

01520  WR ITS ( 8, 5003)  (A VG (KKK) , KKK=1 , NSTEP1 ,5) 

01530  5003  POBBAT  (4X,» AVG», 3X, 20F6. 2) 

01535  WBITE(8,  6667) 

01540  NSS  D=  NS  ED* 1 

01550  to  TO  11 

01560  99  CONTINUE 

01570C  OUTPUT. .. .FINAL  RESULTS 

01580  WRITE (8,5005) 

01590  5005  F03BAT( 1H1 ,//T40,*SUNBATION  OF  SUSPENDED  SEDI RENTS  *) 

01600  WRITE  (8, 5001)  (ADELX  (I)  ,  1=  1,  NST  EP  1,  5) 

01610  DO  412  1=1, N 

01620  WRITE (8,5000)  ADELY (I) , (CSUB ( I, J) , J= 1, NSTEP 1 , 5) 

01630  412  CONTINUE 

01640  WRITE (8,5003)  ( Af GSU B  (KNK) , K KK  =  1 , N STEP  1 , 5) 

01650  WRITE(9, 6668) 

01660  6668  FORRAT (1  HI ) 

01670  DO  411  1=1,21 

01680  AZ=  (1-1) *DELZ 

01690  WR  I TE  (  8,  5002)  A  Z  ,  ( Z  (I  ,  J)  ,  J=1  ,  NSTEP1  ,5) 

01700  5002  FORRAT(4X,F5.2,2X, 20F6. 3) 

01710  411  CONTINUE 

01720  DO  800  11=1,20,4 

01730  YVAL=(IY- 1) *DSLY 

01740  WRITE  (8 ,5555)  TV  AL 

01750  5555  FORRAT ( 1H 1 , ////,T7, *DI STRIBUTI OS  OF  SEDIHENT  IN  BORIZONTAL  *, 
01760*-  ’PLANE  AT  DEPT  H*,  F6  .  2,  *  B  (BG/L)*,//) 

01770  DO  810  IZ= 1 , 21 

01780  IAZ=22-IZ 

01790  AZ=(-DELZ»(IAZ-1)) 

01  BOO  DO  801  IT= 1 , NSTEP1 , INC 

01810  IOUT(IX)  =(CSUB  (I  Y,I  X)  *Z  (I  A  Z,I  X)  *  1000.  ’0 . 5  )/800 

01820  801  CONTINUE 

01830  WB ITS ( 8, 6666)  A Z, (I  OUT ( NKK)  , KKK*1 , NSTEP1 ,5 ) 


01840  6666  FOBEAT  (2X, F6. 2 , 2X, 2016) 

01850  810  CONTINUE 

01860  DO  802  IZ=2,21 

01870  AZ*DELZ*(IZ-1) 

01880  DO  803  IX* 1 , NSTEP1 , ISC 

01890  IOUT(IX)  =  (CSUB  (IT, IX)  *Z  (IZ,IX)  *1  000. +0 . 5 ) /800 

01900  803  CONTI SO  E 

01910  HBITE ( 8, 6666)  A Z, <1 OOT (KKK) , KKK*1 ,NSTEP1 ,5) 

01920  802  CONTINUE 

01930  SHITS  (  8,  6667) 

01940  6667  FOB EAT (//) 

01950  SB  ITS  {  8,  5001)  (ADELX  (KICK)  ,KKR=1  ,  NSTEP1  ,5) 

01960  800  CONTINUE 

01970  STOP 

01930  END 

01990  SUBHOUTINE  A BB A T (BODE ,1 , J, N , 8, S ,D) 

02000C 

020  IOC  CONTESTS  DATA  ABRAY  PROB  SINGLE  TO  DOUBLE  PRECISION  OB  VICE-VERSA 

02020C 

02030C  R0DB=1  -  PBOB  SINGLE  to  dooblb  precision 
02040C  B0DE=2  -  DOUBLE  PRECISION  TO  SINGLE 

02050C  1=  -  BOSS  IN  DATA  BAT8IX 

02060C  J=  -  COLUBNS  IN  DATA  BATRIX 

02070C  N=  -  BOSS  SPECIFIED  IN  DIBSNSION  STATEBENT  FOR  BATRIX  D 

020B0C  fl=  -  COLUBNS  SPECIFIED  IN  DIBENSIC'J  STATEBENT 

02090C  *«**•**«*••«*•*********•************♦****•*«*♦•**•**•*•***••••*.., 

02100  DIBENSION  S(1),D(1) 

02110  NI-N-I 

02120  IF  (BODE-1)  100,  100, 120 

02130  100  IJ  =  1*4 ♦ 1 

02140  NB=N*J»1 

02150  DC*  110  K  =  1,J 

02160  *B=NB-NI 

02170  DO  110  L=1 , I 

02180  IJ=IJ-1 

02190  SB  =NB- 1 

02200  110  D  (NR)=S  (IJ) 

02210  GO  TO  140 

02220  120  IJ=0 

02230  NR  =  0 

02240  DO  130  K=1,J 

02250  DO  125  L=1,I 

02260  IJ=IJ*1 

02270  NB*NR*1 

02280  125  S  (I  J)  =  D  (NR) 

02290  130  SB  =NN  +N I 

02300  140  RETURN 

02310  END 

02320  SUBROUTINE  GRPRD (A, B, R, N, B, L) 

02330C 

02340C  HULTIPLI3S  TSO  RATRICES  TO  FORB  NES  BATRIX 

02350C  *****•*.•••*•••**»•«••*«••*•**.•**•**••••*•*•*•••••••.****.•**•**. 

02360C  A*  -  FIRST  BATRIX 

02370C  B*  -  SECOND  BArTRIX 

02380C  B*  -  OUTPUT  BATRIX 

02390C  N  *  -  ROSS  IN  A 

02400C  fl=  -  COLUBNS  IN  A 

024  IOC  l*  -  COLUBNS  IN  B 

02420C 

02430  DIBENSION  A  ( 1)  ,B  ( 1 )  ,8  (1 ) 

02440  1 8*0 
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02450 
02460 
02470 
02480 
02490 
02500 
02510 
02520 
02530 
02540 
02550 
02560  10 

02570 
02580 
02590 
02600= 
02610C 
02620= 
02630C 
02640= 
02650C 
02660C 
02670C 
02680= 
02690= 
02700= 
02710 
02720 
027  30 
02740= 
02750C 
02760  1 

02770 
02780 
02790 
02800 
02810 
02820 
02830  2 

02840 
02850  3 

02860 
02870= 
02830= 
02890= 

02  900 
02910 
02920C 
02930C 
02  94  0 
02950  4 

02960  5 

02970  6 

02930  7 

02990 
03000 
03010 
03020= 

030  30= 
03040= 
03050 


IK=-8 

DO  10  K=1,L 
IK  =  IK*H 
DO  10  J=1,N 
IB=IR*1 
JI= J-H 
IB*  IK 
R  (I  H)  =0 
DO  10  1=1,8 
JI=JI*N 
IB  =  IB  ♦  1 

R  (IH)  =  R  (IR)  *A(JI)  *B  (IB) 

BETOHS 

BHD 

S0B8DUTIHE  DS2LG  (R  ,  X  ,  8  ,  N  ,EPS  ,1  E  8) 


TO  SOLVE  A  GENERAL  SYSTEM  OF  LI  NBAS  EC'JATIONS 
R  -  DOOBLE  PRECIS  109  8  X  H  RIGHT  HAND  SIDE  HATRIX 
A  -  DOOBLE  PREOISION  8  X  N  COEFF  HATRIX 
8  -  HO  OF  EQUATIONS 

H  -  H0HBE8  OF  RIGHT  HAND  SIDE  VECTORS 
EPS  -  TOLERANCE  POR  TEST 
IER  =  0  -  NO  ERROR 

-1  -  NO  RESULT  BECAUSE  8  LESS  TBAN  1  OR  PIVOT  ELEHENT  =0 
K  -  WARNING  DOE  TO  POSSIBLE  LOSS  OE  SIGNIFICANCE  INDICATED 


DI8ENSI0N  A  (1)  ,R  (1) 

DOOBLE  PRECISION  R,  A, ?IV,TB,TOL, PIVI 
IF (8) 23,23,1 

SEARCH  FOR  GREATEST  ELEHENT  IN  A 
tE3*0 
PIV^O.DO 
88=  8*8 
M8=N*S 
DO  3  L=  1 , 8  8 
TB  =  DA B S  (A  (L)  ) 

IF  (TB-PI7) 3,3,2 

PIV=TB 

I=L 

CON  TINOE 
TOL=  EPS*  PI V 

A  (I)  IS  PIVOT  ELE8ENT 
START  ELIMINATION  LOOP 
LST=  1 

DO  17  K  =  1 , 8 

TEST  ON  SINGULARITY 
IF(PIV)  23,23,4 
IF ( IER)  7,5,7 
IF  (PIV-TOL) 6,6,7 
ISB  =K-  1 

PIVI=1.  DO/A  (I) 

j=(i-i)  /n 

1  =  1 -J*  8-K 
J=J  aI-K 


I*K  IS  ROW  INDEX,  J* K  C0LU8N  INDEX  OF  PIVOT  ELEBENT 

PIVOT  ROW  REDUCTION  AND  ROW  INTERCHANGE  IN  RIGHT  HAND  SIDE  R 

DO  8  L*K,N8,H 


03060 

LL=  L*I 

03070 

TB=PI7I*B (LI) 

03080 

H(LL)  =  R  (L) 

03090 

8 

R(L) =TB 

031 00C 
031 IOC 

IS  ELIHINATION  TEBHINATED 

03120 

IP (K-H)9,18, 13 

03130C 
031 4 OC 

COL0HN  INTERCHANGE  IN  A 

03150 

9 

LEND=LSTtFi-K 

03160 

IP (J) 12,12,10 

03170 

10 

II=J*(1 

03180 

DO  11  L=LST , LEND 

03190 

TB  =  A  (L) 

03200 

LL=L»II 

03210 

A  (L)  =A  (LL) 

03220 

11 

A  (LL)  =  TB 

03230C 

BOS  INTERCHANGE  ADO  PITOT  BOB  BEDDCTICB  IB  A 

03240 

12 

DO  13  L=LST,EH,H 

03250 

LL=L»I 

03260 

TB=  PIVI*  A  (LL) 

03270 

A  (LL)  -A  (L) 

03280 

13 

A  (L)=TB 

03290C 

03300C 

SATE  COLO !1N  INTERCHANGE  INFO 

03310 

A  (L  ST) 

03320C 

03330C 

ELEHENT  REDOCTIOK  AND  NEXT  PITOT  SEARCH 

03340 

PIV=0. DO 

03350 

LST=LST»1 

03360 

J=0 

03370 

DO  16  II=LST,LERD 

0J380 

•PITI=-A  (II) 

03390 

iST=il*n 

03400 

J=J»1 

03410 

DO  15  L=IST,BH,H 

03420 

LL=L-J 

03430 

A  (L)  =A  (L)  tPITI  «A  (LL) 

03440 

TB=  DA 8S  (A  (L)  ) 

03450 

IF(TB-PIV)  15,15,14 

03460 

14 

PIV=TB 

03470 

I  =L 

03480 

15 

CONTINOE 

03490 

DO  16  L  *K,  N  H  ,H 

03500 

LL-  It  J 

03510 

16 

R  (LL)  *R  (LL)  ♦  PI  VI  *R  (  L) 

03520 

17 

LST=  LST  t  fl 

O3530C 

03540C 

END  OP  ELIMINATION  LOOP 

03550C 

BACK  SUBSTITUTION  AND  INTERCHANGE 

03560 

18 

IP  (*-1'  z  3, 22 , 1  9 

03570 

19 

iST=Hn  tn 

03580 

LST=H*1 

03590 

DO  21  1*2, n  - 

03600 

11=  LST-I 

03610 

IST=IST-L3T 

03620 

L=IST-H 

036  30 

L=A  (L)  ♦  0 •  5D 0 

03640 

DO  21  J*II,NH,n 

03650 

TB=R(J) 

03660 

LL=  J 
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03670 

03680 

03690 

03700 

03710 

03720 

03730 

03740C 

03750C 

03760 

03770 

03780 

03785C 

03790C 

03800C 

0381 0C 

03820C 

03870: 

03880C 

03890C 

03900C 

03910= 

03920C 

03930: 

04040C 

04050: 

04060C 

04070: 

04080C 

04090: 

04100C 

04iio: 

04120 

04130C 

04140 

04150C 

04160 

04170 

04180 

041  90* 

04200: 

04210: 

04220 

04230* 

04240* 

04250* 

04260* 

04270 

04280* 

04290* 

04300: 

04310C 

04320 

04330* 

04340* 

04350* 

04360* 

04370* 

04380* 

04390* 

04400 


DD  20  K»IST,HH ,H 
LL=  LL* 1 

20  TB  =  TB-A  (K)  *B  (LL) 

K=  J*L 

S(J)  =R  (N) 

21  S(K)=TB 

22  RETURN 

ERROR  RETORN 

23  IER=-1 
BETUBN 
END 

IHSL  ROOTIRE  NAHE  -  flERF=ERF 


PURPOSE 

-  EVALUATE  THE  ERROR  FUNCTION 

USAGE 

-  RESULT  =  ERF(Y) 

A  RGU  RENTS 

u. 

OS 

-  INPUT  ARGUMENT  OP  THE  ERROR  FUNCTION. 

-  OUTPUT  VALUE  OF  THE  ERROR  FUNCTION. 

COPYRIGHT 

-  1978  BY  I  SSL ,  INC.  ALL  BIGHTS  RESERVED. 

WAB6ANTY 

-  IHSL  WARRANTS  ONIY  THAT  I HS L  TESTING  HAS 
APPLIED  TO  THIS  CODE.  NO  OTHER  WARRANTY 
EXPRESSED  OR  IMFIfSD,  TS  APPLICABLE. 

REAL  FUNCTION  ERF(T) 


SPECIFICATIONS  FOR  ABODE ENTS 


BEAL 

INTEGER 
DI SESSION 
REAL 


DATA 


DATA 


DATA 


DATA 


T 


SPECIFICATIONS  FOR  LOCAL  VARIABL 

ISB,I 

P(5),C(3),P1(S),C1(7)  ,?2(S)  ,2  2(4) 

P,Q,P1 ,01 ,P2 ,22 ,XMIN ,XLARGE,SSQPI,X, 
RES,XSO,XNUM,XDES,XI 

COEFFICIENTS  FOR  0.3  . LE.  Y  .LT. 
.  477 

P(1)/-. 4 44 226 47396874/, 

P (2) /10. 731707253648/, 

P (3) /15. 915606197771/, 

P  (4) /374. 8 162400 1284/, 

P (5) /2.5612422994323E-02/ 

Q(1)/17. 903143558843/, 

0 (2) /124 .82P92031531/, 

Q(3j/332.  17224470532/ 

COEFFICIENTS  FOB  .477  .LE.  Y 
.LE.  4.0 

PI (1J/7.211758250R831/, 

PI (2) /43. 162227222357/, 

PI (3)/152. 98928504694/, 

PI (4) /339, 32031673434/, 

PI (5)/451.918953711R7/, 

PI (6J/300. 45926102016/, 

PI (7) /- 1 . 3636485738272E-07/, 

PI  (8)/.  56419551747897/ 

01 (1)/77. 00U152935229/, 


04410* 

Q1  (2J/277. 58544474399/, 

04420* 

01 (3)/638. 98026446563/, 

04430* 

01  (4) /931.  3540948506  1/, 

04440* 

01  (5)/790. 95092532790/, 

04450* 

01  (6J/300. 45926095698/, 

04460* 

Q1  (7) /1 2. 782727319629/ 

044700 

COEFFICIENTS  FOR  4.0  .LT. 

04480 

DATA 

P2  (1)/-. 22695659353969/, 

04490* 

P2  (2)/- 4. 947 309 10 62 3 252- 02/, 

04500* 

P2  (3)/-2. 99610707703542-03/, 

04510* 

P2  (4) /-2.  2319  2459734  182-02/, 

04520* 

P2 (5)/- 2. 7866130860 9 65E-01/ 

04530 

DATA 

Q2 (1) /I. 05 1675 107067 9/, 

04540* 

Q2 (2)/. 19136892610783/, 

04550* 

Q2  (3) /I. 06269  230  528  4  72-02/, 

04560* 

Q2 (4) /I . 9873320131714/ 

045700 

CONSTANTS 

04580 

DATA 

XHIN/1 .02-8/, X  LARS  2/5. 63  75 EO/ 

04590 

DATA 

SSQPI/. 56418953354776/ 

04600C 

FIRST  EXECUTABLE  STATEHENT 

04610 

X  =  T 

04620 

ISS  =  1 

04630 

IF  (X.G2.0.0E0)  GO 

TO  5 

04640 

IS»  =  -1 

04650 

X  =  -X 

04660 

5 

IF  (X.LT.. 47720)  GO 

TO  10 

04670 

IF  (X.  12.  4.0  SO)  GO 

TO  25 

04680 

IF  (X.LT.XLARG2)  GO 

TO  35 

04690 

H2S  =  1 .  20 

04700 

GO  TO  50 

047100 

ABSfY)  .IT.  .477,  ETALUATE 

04720C 

APPROXIMATION  FOR  ERF 

04730 

10 

IF  (X.  IT.  X  HI  N)  GO  TO  20 

04740 

xso  =  X*X 

04750 

XNUH  *  P  (5 ) 

04760 

DO  15  1=1,4 

04770 

XNUH  =  XNUM»XSQ* 

P(I) 

04780 

15 

CONTINUE 

04790 

XDEN  =  (  (Q  (1)  *XSQ)  • 

XSQ*Q  (2J )  *Xr,0*Q(3) 

04800 

RES  =  X*XN  UH/XDEN 

04810 

GO  TO  50 

04820 

20 

RES  =  X*?  (4)  /2  (3) 

04830 

GO  TO  50 

04840C 

.477  . IE.  ABS  (Y)  .12.  4.0 

04850C 

EVALUATE  A PP30 X  IE A TION  FOR 

04860 

25 

xsg  *  X*I 

, 

04  870 

XNUH  =  PI  (7 )  * X  * P 1  (9) 

04880 

XDSN  =  X*Q  1  (7) 

04890 

DO  30  1=1,6 

04900 

XNUH  =  XNUH  * X* P 1  (I) 

04910 

XDEN  =  XDEN*X*Q1 (I) 

04920 

30 

CONTINUE 

04930 

RES  =  XRUH/XDEN 

04940 

GO  TO  4  5 

04950C 

4.0  .LT.  ABS(Y) ,  EVALUATE 

04960C 

APPROXI HATION  FOB  ERF 

04970 

35 

XSQ  =  X*X 

04980 

XI  =  1.020/XSO 

04  990 

XNUH  =  P2 (4) *XI*P2 (5) 

05000 

XDEN  =  XI*Q2(4) 

05010 

DO  40  1=1,3 

189 


05020 

XHUH  *  XNUH*XI»P2 (I) 

05030 

XDEN  »  XDEN*XI*Q2 (I) 

050  40 

40 

CON TIH OE 

05050 

RES  *  (SSQPI+XI*XN3M/XDEN) /X 

05060 

45 

BES  *  BSS*EXP (- XSQ) 

05070 

RES  *  1.0E0-RES 

05080 

50 

IF  (ISB.E3 1)  BES  =  -BES 

05090 

ERF  *  BES 

05100 

RETUB  11 

05110 

BND 

I 


APPENDIX  F 

HORIZONTAL  SEDIMENT  DISTRIBUTIONS 


AT  DEPTH  OF  1.2  m 
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-50.00 

0 

c 

0 

0 

0 

0 

0  0 

0 

0 

0 

-47.50 

0 

0 

0 

p 

0 

0 

0  0 

0 

0 

0 

-45. 00 

0 

c 

0 

0 

0 

0 

0  p 

0 

0 

0 

-42.50 

0 

0 

0 

0 

0 

0 

0  ') 

0 

0 

0 

-40.00 

C 

c 

0 

0 

0 

0 

0  0 

0 

0 

0 

-37.50 

0 

0 

0 

0 

0 

0 

0  0 

0 

0 

0 

-35.00 

0 

c 

0 

0 

0 

0 

0  0 

0 

0 

0 

-32.50 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-30. 00 

c 

c 

0 

0 

0 

0 

0  0 

0 

0 

0 

-27.50 

0 

0 

0 

0 

0 

0 

0  0 

0 

0 

0 

-25.00 

0 

c 

0 

0 

0 

0 

0  0 

0 

0 

c 

-22.50 

0 

0 

0 

0 

P 

0 

0  o 

0 

0 

o 

-20.00 

c 

c 

0 

0 

0 

0 

C  L 

0 

0 

0 

-17.50 

0 

0 

0 

0 

0 

0 

P  5 

0 

c 

0 

-15. 00 

c 

c 

0 

0 

1 

2 

3  '4 

5 

6 

6 

-12.50 

1  25 

37 

33 

33 

37 

31 

31  3  • 

30 

7  ) 

77 

-  10. 00 

125 

79 

65 

65 

63 

60 

55  5? 

55 

53 

57 

-7.50 

125 

77 

67 

66 

64 

6  3 

67  61 

60 

5w 

58 

-5.  00 

1  25 

77 

67 

66 

65 

6  3 

62  6’ 

60 

5  8 

5  3 

-2.50 

1  25 

7  7 

67 

66 

65 

63 

62  61 

60 

5  > 

58 

0.00 

125 

75 

67 

66 

65 

63 

6?  6 1 

60 

5  ) 

58 

2.50 

125 

77 

67 

66 

65 

63 

62  61 

66 

5  * 

58 

5.  00 

125 

77 

67 

66 
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